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ABSTRACT 

The  response  of  a  hollow  circular  cylindrical  shell  of  arbitrary  thickness,  in 
either  an  elastic  or  a  viscoelastic  medium,  to  transient  dilatational  and  shear 
waves  (and  their  s  >erposition)  is  presented.  The  solution  is  valid  within  the 
scope  of  the  linear  theory  of  elasticity  or  viscoelasticity.  The  technique  for 
obtaining  the  solution  relies  upon  1)  the  construction  of  a  train  of  incident  pulses 
from  steady- state  components,  where  each  pulse  represents  the  time  history  of 
the  transient  stress  in  the  incident  '*>ave,  and  Z)  the  existence  of  a  physical 
mechanism  Inat,  between  pulses,  restores  the  disturbed  particles  of  the  cylinder 
and  the  surrounding  medium  to  an  unstrained  state  of  rest. 

7b*  jn  the  cylinder  response  of  the  following  factors  is  discussed; 

liner  thickness,  cylinder-medium  impedance  mismatch,  viscoelasticity  in  the 
medium,  and  incident  wave  form  4  step  pulse,  rectangular,  triangular,  linear 
rise- exponential  decay), 
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The  study  reported  herein  is  concerned  with  the  dynamic  re* 
spouse  of  deep- bur led  protective  installations  subjected  to 
stress  waves  resulting  from  nuclear  explosions.  As  pointed  out 
elsewhere,  for  example  in  Ref.  1,  the  actual  problem  is  exceed¬ 
ingly  complex.  However,  certain  justifiable  idealizations  can 
be  made,  that  lead  to  tractable  formulations  of  the  subject 
problem,  the  solution  of  which  yields  information  on  the  princi¬ 
pal  effects  in  the  actual  problem.  For  example,  the  structure 
considered  is  a  long  thick  isotropic  elastic  cylinder.  The  sur¬ 
rounding  medium  is  regarded  in  two  distinct  phases  of  the  study 
as  either  elastic  or  viscoelastic.  It  is  assumed  that  the 
pressure  front  generated  by  the  nuclear  burst  travels  on  the 
surface  with  super  seismic  velocity  and  thus  can  be  considered 
to  transmit  plane  seismic  waves  Into  the  ground  {Ref.  1).  Be¬ 
cause  :  of  the  depth  at  which  the  structure  Is  burled,  the  influ¬ 
ence  of  Rayleigh  waves  Is  neglected.  For  the  same  reason,  it 
is  assumed  that  no  appreciable  influence  is  exerted  on  the  peak 
response  by  reflections  from  the  surface.  Gravity  stresses  in 
the  cylinder  are  not  considered. 


In  view  of  the  foregoing  considerations,  the  idealized  prob¬ 
lem,  which  is  the  subject  of  the  analytical  and  computational 
investigation  described  in  this  report,  is  that  of  determining 
the  stresses  and  displacements  in  an  infinitely  long  hollow 
cylinder  of  arbitrary  thickness  embedded  in  an  infinite  medium. 
The  cylinder  is  enveloped  by  plane  stress  waves  progressing  in 
a  direction  perpendicular  to  its  axis  {see  Fig.  1) .  The  problem 
is  thus  one  of  plane  strain.  The  cylinder  is  elastic,  homoge¬ 
neous,  and  Isotropic.  The  surrounding  medium  is  homogeneous  and 
isotropic,  and  In  two  separate  and  distinct  phases  of  the  analy¬ 
sis  Is  considered  to  be  either  elastic  or  viscoelastic.  The 
plane  waves  impinging  on  the  cylinder  are,  in  succession,  di La¬ 
ta  tional  and  distortions  (Fig.  1).  An  analysis  is  developed 
to  account  for  an  arbitrary  pressure-time  history  of  the  incom¬ 
ing  waves. 


Computations  have  been  carried  out  for  the  following  wave 
forms:  1)  a  step  pulse;  2)  a  rectangular  wave;  3)  a  triangular 
wave;  and  4)  a  linear  rise- exponential  decay  wave  form. 
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Two  limiting  cases  of  the  present  problem  have  been  studied 
by  a  number  of  investigators.  Hie  diffraction  of  a  stress  wave 
by  a  cylindrical  cavity  in  an  infinite  medium  is  examined  in 
Refs.  2,  3,  4,  and  5.  Using  thin  shell  theory  to  describe  the 
behavior  of  the  liner.  Ref.  6  trea;.3  the  problem  of  the  diffrac¬ 
tion  of  a  stress  wave  by  a  thin  cylindrical  liner  in  an  elastic 
medium. 

In  the  foregoing  references,  transform  or  related  integral 
techniques  are  employed  to  obtain  solutions.  Such  techniques, 
however,  if  applied  to  the  present  problem,  where  the  cylinder 
must  be  analyzed  within  the  scope  of  the  theory  of  elasticity 
(since  it  is  considered  to  be  of  arbitrary  thickness),  would 
lead  to  formidable  difficulties  in  the  inversion  process. 

The  method  utilized  in  the  current  analysis  circumvents 
these  difficulties.  It  consists  of  devising  a  solution  to  the 
transient  problem  from  a  superposition  of  appropriate  steady- 
state  solutions  to  the  wave  equation.  Its  success  stems  from: 

1)  the  capability  of  constructing  a  train  of  pulses  from  steady- 
state  sinusoidal  components  (Fig.  2a)  where  each  pulse  represents 
with  sufficient  accuracy  the  time  history  of  the  transient  stress 
in  the  incident  wave;  2)  the  existence  of  a  mechanism  which, 
between  pulses,  restores  the  disturbed  particles  of  the  cylinder 
aud  the  surrounding  medium  to  an  unstrained  state  of  rest. 
Physically,  in  the  present  problem,  this  mechanism  is  provided 
by  the  continuous  rauiation  of  energy  from  the  excited  cylinder 
through  the  surrounding  medl”*"  jutward  to  infinity.  Mathemati¬ 
cally,  it  is  implied  by  solutions  of  the  field  equations  that 
enforce  continuity  of  stresses  and  displacements  between  the 
cylinder  and  the  embedding  medium,  and  which  at  infinity  repre¬ 
sent  cutgoing,  decaying  waves.  Thus,  if  time  is  measured  from 
the  moment  of  arrival  of  a  pulse  at  the  cylinder,  the  response 
will  oe  identical  to  that  excited  by  a  transient  loading  having 
the  time  history  of  a  single  pulse.  The  validity  of  this  tech¬ 
nique  is  amply  demonstrated  in  the  present  report. 

The  use  of  steady-state  solutions  to  obtain  a  transient 
response  has  also  been  suggested  by  Ref.  8  for  systems  in  which 
damping  provides  the  mechanism  for  bringing  the  system,  to  rest. 
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As  will  be  seen,  the  present  approach  offers  greater  opportunity 
to  gain  physical  insight  into  the  problem  than  do  previous  ap- 

pi. uac’una .  An  auultiuiial  point.  u£  uepai  Luie  in  Lne  present 

approach  concerns  the  number  of  circumferential  modes  chosen  to 
describe  the  deformation.  *  11  the  efore“msntioned  approaches 
constrain  the  deformation  by  selecting  in  advance  a  fixed  number 
of  modes  to  apply  to  all  material  points  and  at  all  times  of  the 
response.  In  the  pfetent  analysis,  the  accuracy  with  which  the 
solution  is  required  to  converge  determines  the  number  of  modes 
selected.  This  number  is  permitted  to  vary  with  both  the 
material  point  being  considered  and  the  time  at  which  its  re- 
sponse  is  calculated.  Thus,  within  the  limitations  of  a  Fourier 
series  representation  at  the  initial  discontinuity,  the  present 
method  should,  in  general,  be  more  accurate.  This  will  be  shown 
to  be  especially  the  case  during  the  early  part  of  the  response. 

The  dynamic  response  of  the  cylinder  is  determined  by  ob“ 
twining  separately  and  subsequently  combining  several  distinct 
component  parts,  each  e£  which  is  a  solution  of  the  wave  equation 
in  either  the  cylinder  or  the  surrounding  infinite  medium.  These 
are  the  stresses  and  displacements  associated  with:  1)  the  inci- 
dent  wave;  2)  the  transmitted  and  reflected  waves  in  the  cylinder 
and  3)  the  scattered  waves  in  the  surrounding  medium.  The  fore- 
going  components  must  be  combined  so  that  the  inner  surface  of 
the  cylinder  is  traction  free,  and  at  the  outer  surface  the  normal 
stresses  in  the  radial  direction,  the  shear  stresses  in  the  cir¬ 
cumferential  direction,  and  the  displacements  in  the  cylinder  are 
continuous  with  the  corresponding  quantities  in  the  medium. 
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I.  Stresses  and  Displacements  Applied  to  the  Cylinder  by  an 
Incident  Dilatational  Plane  Wave  in  an  Elastic  Medium 

The  Cartesian  coordinate  system  in  which  the  incident  plane 
dilatational  wave  is  described  is  illustrated  in  Fig.  3a.  The 
disturbance  is  assumed  to  be  traveling  in  the  negative  x- direct ion. 
The  equation  governing  the  propagation  of  the  wave  (which  depends 
only  on  the  space  coordinate  x  and  time  coordinate  t)  is 
given  by  Ref.  9  as 


4X 


(3.1.1) 


where  4>  is  a  displacement  potential  defined  as  follows: 


x  dx 


uy  ■  0 
u*  ■ 0  > 


(3.1.2) 


and  the  pertinent  stresses  are  given  by 


o  -  (Tv  +  2u  )  — ~  , 
xx  v  m  Har  ^2  * 


(3.1.3) 
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0  m  A  , 

yy  m  av2 


(3.1.4) 


The  quantities  Ujj,  uy,  u2,  o^,  and  o  are 
and  stresses  in  Cartesian  coordinates,  Ajq  anr 
Lame  constants  of  the  medium,  c^  is  the  spe 
of  dilatational  waves  in  the  elastic  medium  anc 


raents 

ne 

;ation 

fcy 
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cdm 


+  2n. 
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m 


(3.1.5) 


»nij  v  1  q  t-he  mass  density, 
‘m 


A  steady  state  solution  of  (3.1.1)  is 


$  ■  Y  A  e 
L  p 

P-1 


i[kp3(x-b)  +  aipt] 


(3.1.6a) 


where  the  frequency  cd^  is  given  by 


and 


El 

to  ■  „ 
P  T 


.a. 


p3  c 


dm 


(3.1.6b) 

(3.1.6c) 


Also,  T  is  an  arbitrary  time  interval  and  the  Ap  are  complex 

constants  of  superposition  that  will  be  expressed  in  terms  of  the 
Fourier  coefficients  defining  the  incoming  wave  form. 

For  our  purpose,  we  represent  the  incident  stress  wave  o  , 
(see  Fig.  2)  in  the  form 


a  ■ 

xx 


l 

p*»l 


a  sin  u> 
P  P 


(*rf> ♦  * 


*  * (3.1.7a) 


or  at  x  ■  b  as 


xx  L,  p 
P-1 


Za  sin  o>  t  , 
P  P 


(3.1.7b) 


where  the  ap  are  the  coefficients  defining  the  wave  form  over 
the  arbitrary  half  period  T. 
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3  *  T 

P  T 


a  sin  cu  t  dt 

XX  p 


(3.1.7c) 


pare  of 


An  expression  for  identical  with  (3.1.7a)  is  the  real 


a  m 

xx 


1  <•  UP> a 

p-1 


+copt] 


(3.1.8) 


From  (3.1.3)  and  (3.1.6a), 
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Con^aring  (3.1.8)  and  (3.1.9),  we  obtain 


(3.1.9) 
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(3.1.10) 


Now,  from  (3.1.2),  (3.1.6a),  and  (3.1.10)  we  can  write 
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tlkp3(*-b) 


v1 


or,  noting  (3.1.5)  and  (3.1.6c), 


(3.1.11a) 
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(3.1.11b) 
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Since  the  cylinder  boundaries  are  naturally  defined  by  polar 
coordinates  it  will  be  more  convenient  to  write  the  stresses  and 
displacements  associated  with  the  incident  wave  in  oolar  coordi¬ 
nates.  Then 


°rr  "  *xx(cosi£e  +  €  sin2°)  >  (3.1.12a) 


<4  "  "  CTxx  (H"*)  aln  20  '  (3.1.12b) 

where 

X 

e  "  r •  (3.1.12c) 

m  m 

Substituting  (3.1.8)  into  (3.1.12)  and  setting  x  -  b  cos  9, 
we  obtain  the  following  expressions  for  the  stresses  applied  to 
the  boundary  by  the  incident  wave: 

00 

jn  icu it 

rr  "  ■  2j  °x,  p  e  *  (3.1. 13a) 

P-1 


B  ^  io>  t 

are  "  ■  1  are,p  e  '  (3.1.13b) 

P-1 

where 

2  2  iq-(cose-l) 

ax,p  “  ("iap)(cos  ©  +  €  sin  ©)  e  ,  (3.1.13c) 


,  v/l  -  ex  iq3(cos  ©-1) 

°r©,p"  ('iapH  2  )8ln  20  e  >  (3.1.13d) 
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with 


q3  •  bkp3  "  Krjp)  *  (3.1. 13®) 

dm 

In  anticipation  of  matching,  later  in  tha  analysis,  the 
tractions  applied  by  the  incident  train  of  pulses  with  the 
cylinder  stresses  and  stresses  associated  with  scattered  waves 
in  the  medium,  we  represent  the  stresses  given  by  (3.1.13), 
and  acting  at  the  cylinder  boundary,  as  the  real  part  of 
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(3.1.14a) 
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(3.1.14b) 


where  the  complex  Fourier  coefficients  Apn  and  are 

obtained  by  comparing  (3.1.14)  with  (3.1.13).  They  are 
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n  £  1 


A  n  “ 

pO 


1 

TT 


2  2  iq3(cos  0-1) 

(cos‘0  4-  e  sin*0)  e  J  d 9  , 


(3.1.15b) 
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(3.1.15c) 
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These  expressions,  (3.1. IS),  are  representable  in  terms  of 
Bessel  functions  of  a  complex  argument  of  the  first  kind  and 
integral  order  n  —  as  shown  in  Appendix  B. 

Conditions  of  continuity  between  the  cylinder  and  surround¬ 
ing  medium  require  that  we  match  displacements  as  well  as  appro¬ 
priate  stress  components.  The  radial  and  tangential  displacements 
in  polar  coordinates  are  related  to  the  Cartesian  coordinate 
displacement  U*  as  follows: 


ur  -  ux  cos  6  ,  u0  - 


u  sin  9 
x 


Substituting  (3.1.11b)  into  (3.1.16)  and  setting 
x  ■  b  cos  9  leads  to 
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(3.1.16) 


(3.1.17a) 


(3.1.17b) 


(3.1.17c) 


(3.1.l7d) 


or  alternatively. 
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As  with  Che  stresses,  we  expand  the  displacements  at  the 
boundary  in  the  form 
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where  the  complex  Fourier  coefficients  and  Dpn  are 

obtained  by  comparing  (3.1.18)  with  (3.1.17).  Thus, 


pn  tt 


iq3(cos  9-1) 

e  cos  9  cos  n6  d9  , 

0  n  >  1 


(3.1.t9a) 


°P0- 


1 

T 


1^3 (cos  9-1) 


cos  9  d9  , 


(3.1.19b) 


Dpn  " 


l 

IT 


iq3(cos  9-1) 


sin  9  sin  n9  d9 


(3.1.19c) 


As  in  (3.1.15),  the  complex  Fourier  coefficients  defined 
by  (3.1.19)  are  shown  in  Appendix  B  to  be  representable  in  terms 
of  Bessel  functions  of  a  complex  argument. 
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Displacements  Uj.  and  u q,  which  satisfy  (3.2.1)  through 
(3.2.3)  Identically,  and  are  associated  with  dilataCional  waves, 
are  given  by 


$ 


d 

U2 


JL 

_  * 

I 


(3, 2,4) 


in  which  $c,  the  displacement  potential,  is  governed  by  the 
wave  equation 


2  1 

rte  -  ,  — f-  -  o  . 

c  Cj  dt 

o 


(3.2.5a) 


where 


„2 .  aL  +  i  a.  +  jl  af¬ 
ar2  *  *  r2  d>2  ' 


<3. 2. 5b) 


and  the  speed  of  propagation  of  dilatational  waves  in  the 

cylinder,  is  given  by 


/THE 

V  v  * 


(3.2.5c) 


where  y  is  the  mass  density  of  the  cylinder. 

When  Eqe.  (3.2.4)  are  introduced  into  the  strain  displace* 
went  relations  (3.2.1),  and  these  in  turn  into  the  stress  strain 
rwationa  (3.2.2),  we  obtain  the  stresses  associated  with  dila- 
tat tonal  waves,  in  terms  of  the  displacement  potential  as 

follows 


rr 


ts 


•  XV V  +  2;i 


bh 


Sr 


°ee 


o  ,  i  d2<*>c  i 

*7‘4>o  + 


-  2|i 
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h  ■ 


(3.2.6) 
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A  solution  of  (3.2.5a)  may  be  written  as 


go  (X 

f  r  1  >  t 

®c-  11  1  K.pnVV0  +Z2,pnVV>jc°sne}e  P 

o*«  1  n— fl 
* 

(3.2.7) 

00 

+  y  Ap^r  cos  9  , 

p-1 


where:  1)  k  ^  the  series  )  Aplr  cos  9  cons*-sts 

p-1 

of  a  superposition  of  solutions  to  (3.2.5)  that  correspond  to  rigid 

body  displacement;  3)  the  Jn(«)  and  Yn(«)  are  Bessel  functions 

of  the  first  and  second  kind,  respectively,  of  integral  order  n; 

and  4)  the  Z,  and  Z„  are  complex  constants  of  super - 
Ijpn  Z,pn 

position  to  be  determined  from  the  boundary  conditions. 

A  bar  is  placed  over  the  integer  n  in  certain  subsequent 
terms  for  the  following  reason:  In  the  analysis  to  follow,  the 
potential  functions  will  be  differentiated  with  respect  to  the 
circumferential  coordinate,  9,  and  also  the  radial  coordinate, 
r.  Differentiation  with  respect  to  9  introduces  the  integer  n 
as  a  factor  directly  into  the  expressions  for  the  stresses.  Dif¬ 
ferentiation  with  respect  to  r,  and  subsequent  use  of  the  Bessel 
function  recursion  relations,  also  introduces  the  integer  n  into 
the  stress  and  displacement  expressions.  It  will  be  convenient 
to  identify  the  origin  of  the  integer  n  because  the  solutions 
for  the  cylinder  and  for  the  scattered  response  in  the  medium, 
appropriate  for  incoming  shear  waves,  can  be  obtained  from  the 
corresponding  solutions  for  incoming  diiatational  waves  by: 

1)  interchanging  sin  nd  and  cos  nS;  and  2)  setting  H  »  -  n. 

In  the  case  of  the  diiatational  wave,  we  simply  set  n  -  n. 

It  will  be  noted  that  $c  is  constructed  by  a  superposition 
of  steady  state  solutions  of  the  form 
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and  a  solution  describing  rigid  body  displacements.  Such  solutions 
have  been  used  to  solve  the  steady  state  problems  reported  in 
Refs.  11  and  12.  The  use  of  steady  state  solutions  to  construct 
the  solution  to  a  transient  problem  has  already  been  discussed. 

tt  is  shown  in  Ref.  13  that  the  derivatives  of  the  Bessel 
functions  may  be  expressed  as 


i 


(z) 


“  J  (z) 
z  n  w 


Jn+l(z)  ' 


(3.2.9a) 


JnU) 


n(p,  :.,1)  .  1 


Vz)  +  z  -Wz>  ' 


(3.2.9b) 


and  similar  relations  for  the  derivatives  of  Yn(z)  • 

As  mentioned  in  the  text  preceding  Eq.  (3.2.8),  and  as 
can  be  seen  in  Eqs.  (3.2.9),  the  recursion  relations  for  the 
derivatives  of  the  Bessel  functions  introduce  the  integer  n 
into  the  stress  and  displacement  expressions.  Introducing 
Eq.  (3.2.7)  into  Eqs.  (3.2.4)  and  (3.2.6),  making  use  of  rela¬ 
tions  (3.2.9a  and  b),  and  defining  a  nondime.-sional  radial 
coordinate  as 


P  -  £  ,  (3.2.9c) 

and  a  nondimens ional  wave  number  as 

qx  -  bkpl  ,  (3.2.9d) 

wa  can  write  the  expressions  for  the  displacements  and  stresses 
as 


00  00 


Ur  “  b  1  b  ql{Zl,pn(q"p  Jn(qlp)  "  Jn-41(qlp)) 


p-1  n«0 


+  Z2,pn(^  W>  ‘  Wqlp>)}c03  1,6  e 


iui  t 

p 


A  ,cos 
Pi 
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(3.2.10a) 
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a>  cl 
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00 


CO 


'«  '  72  l  l  ;%{Zl.pnivJn+l(V>  *  <n  •  L>Jn<V>j 


(—>  l—J  - 

o 

p—  X  II— JL 


(2 .2. lOe) 


+  Z2,pn[vWV)  '  (n  ‘  1>W>]}  sln  "e 


io>  t 

e  P  , 


where  the  superscript  d  has  been  used  to  designate  the  con¬ 
tribution  from  the  dilatational  solution.  The  contribution  from 
the  distortional,  or  shear  wave  solution,  will  now  be  discussed. 

A  representation  of  the  displacements  u*.  and  u0  associa¬ 
ted  with  shear  waves  is 


r  r  d@ 


__l£ 
Sr  ' 


(3.2.11) 


in  which  the  displacement  potential,  1  ,  is  governed  by  the 
wave  equation  c 


_2  1_  0  *c.  A 

*c  ct  dt2 


(3.2.12a) 


where  ct,  the  speed  of  propagation  of  shear  waves  in  the 
cylinder,  is  given  by 


ft 

\  "  A  * 


(3.2.12b) 


When  Eqs.  (3.2.11)  are  Introduced  into  the  strain  displacement 
relations  (3.2.1)  and  these  in  turn  into  the  stress  strain 
relations  (3.2.2),  we  obtain  the  stresses  associated  with  shear 
waves  in  terms  of  the  displacement  potential  1?  : 


■N  ,  OI 

iL  1 _ s 

dr  r  dS 


v  L  df 

_  r  i_  I _ s 

r  dr  r  dr 


'  crr  ' 


| 


(3.2.13) 


r2  de2 
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A  solution  of  (3.2.12a)  may  be  written  as 


I  {  I  [Z3,pnW>  +  Z4,pnWr> 

p-1  n-1 


A  « 

a in  n©  l  e  ^ 


00 


+  £  Bplr  sin  e  , 

p-1 


(3.2.14) 


where:  1)  kp2  -  cop/ct;  2)  the  series  ^  ®plr  3*-n  0  con- 

p-1 

sists  of  a  superposition  of  solutions  to  (3.2.12)  that  corre¬ 
sponds  to  rigid  body  displacements;  and  3)  the  Z,  and 

Z4  pn  are  complex  constants  to  be  determined  from  the  boundary 
conditions . 

As  in  the  dliatational  solution.  ¥  is  constructed  from 

c 

a  superposition  of  steady-state  solutions  of  the  form 


% 


l  [ 


n« 


Z3,pnWr) 


+  Z4,p„W' 


sin  n0  e 


io>  t 
P 


(3.2.15) 


and  a  solution  corresponding  to  rigid  body  displacements. 

When  Eq.  (3.2.14)  is  substituted  into  Eqs.  (3.2.11)  and 
(3.2.13),  and  the  relations  in  (3.2.9a  and  b)  and  the  quantity 
defined  in  (3.2.9c)  are  introduced  along  with 


q2  -  bkp2  ,  (3.2.16) 

we  have  finally  the  following  expressions  for  the  displacements 
and  stresses  associated  with  the  shear  wave  solution: 
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I  I  ?{Z3,pnVV>  +Z*,p»VV>} 

3-1  n-0  1 


leu  t 

cos  nd  e  P 


I  Bpl  cos  0  > 


P-1 


(3.2.17a) 


00  00 


b  ^  ^  q2{Z3,pn 

P-1  n-1 


q2P  jn(q2p)  “  Jn+l(q2p> 


+  Z^,pn  q^p  Yn^q2p)  *  Yn+l^q2p^ 
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io)  t 
sin  nS  e  p 


(3.2.17b) 


I  Bpl  sin  6  , 


P-1 


00  00 


°’r  '  S  \  \  {Z3,pn[<n  -  1>W>  -  VWV>] 

p— 1  n— 0  r 


(3.2.17c) 
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ico  t 


(n  -  l)Yn(q2p)  -  q2pYn+i(q2p) jjeos  nS  e  p 


a3  m  -  rrs 

qQ0  arr  * 


(3.2.17d) 
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+  r^r  J. .  (q„p) 

HoH  irrA  *■ 
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4,pn 


/  2n(n  -  1)  .\„  , 

(“72  "  WV' 


(3.2. 17e) 


(q2P) 


+  v  Wv> 


} 


ia:  t 
sin  nS  e  ^ 


3.  Solution  of  the  Wave  Equation  for  Scattered  Waves  in  the 
Surrounding  Medium 

As  in  the  cylinder,  the  solution  in  the  surrounding  infinite 
medium  may  be  constructed  from  two  displacement  potentials,  one 
dilatational,  the  other  distortional.  However,  in  addition  it  is 
required  that  at  infinity  these  solutions  represent  outgoing, 
decaying  waves. 

After  04  is  replaced  by  c^  (the  dilatational  wave 
velocity  in  the  medium) ,  the  dilatational  potential,  namely  a 
solution  of  (3.2.5a)  that  meets  the  requirements  at  infinity, 
may  be  taken  as 


00  00  ,  ^  00 

(o\  10)  t 

\m  1  1  Z5,pnHn  <kP3r) CO8n0  6  P  +  I  Cplr  cos  0  ,  (3.3.1) 

p-l  n-0  p«l 

where  l)  k  i  2)  H<2)  -  J  (*)  -  «.(») 

^  dm 

are  Hankel  functions  of  the  second  kind;  3)  the  Zc  „  are 

5,pn 

complex  constants  of  superposition  to  be  determined  from  the 
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00 


boundary  conditions;  and  4)  the  series  Y  ^pl*  r  cos  9  con- 

p-1 

sists  of  solutions  associated  with  rigid  body  displacements. 

Substituting  Eq.  (3.3.1)  into  Eqs.  (3.2.4)  and  (3.2.6)  and 
introducing  Eqs.  (3.2.9a  and  b)  and  the  definition  (3.2.9c)  along 
with  the  definition  q^  «  we  obtain  the  following  expres¬ 

sions  for  displacements  and  stresses: 


“f  -  i  I  l  Kpn[?  W> 

p-1  n-0 


q3Jn+l 


(q3p) 


1Z5,pn[?  Yn(q3p>  '  q3Yn+l<q3p) ]}cos  n0 


(3.3.2a) 


l  cpi  c°8  8  - 


oo  oo 

’  b  I  I  {z5,pn[p  Jn(q3p).  "  iZ5,pn[p  Yn(q3p)]}sln  n0 


ioi  t 
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(3.3.2b) 
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(3.3.2c) 
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00 
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where,  as  already  cited, 

<13  -  t>“p3  •  (3.3.2e) 

After  ct  is  replaced  by  ctm  (the  distortional  wave 
velocity  in  the  medium) ,  the  distortional  displacement  potential 
for  the  medium,  namely,  a  solution  of  (3.2.12a)  that  nreets  the 
requirements  at  infinity,  may  be  taken  as 


~  (o\  io>t  _ 

*■*  1  l  Z6,pnHn  Vr)Sl"ne  e  +  l  Dplr  alr‘  6  >  <3-3-3> 

p-t  n-1  p»l 

where:  1)  the  series  £  Dpir  sin  9  consists 

of  solutions  associated  with  rigid  body  displacements;  and 
3)  the  are  complex  constants  of  superposition  to  be 

determined  from  the  boundary  conditions. 

The  displacements  and  stresses  associated  with  ¥  are 

m 

obtained  by  substituting  (3.3.3)  into  equations  of  the  form 
(3.2.11)  and  (3.2.13),  and  introducing  (3.2.9a-c).  Thus 


ur"-b  l  l  ?{Z6.Pn[W>] 


p-1  n«0 


(3.3.4a) 
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(3.3.4b) 
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(3.3.4c) 
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ioo  t 
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where 


**4  "  bkp4  ‘ 


(3.3. 4e) 
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4.  Correspondence  Between  Elastic  and  Viscoelastic  Field  Equations 
Viscoelastic  problems  are  very  often  formulated  on  the  basis 

.  ^  -  1 i  <  •  _  i  »  i  .*  i  ^  v  • 

wj.  «u  auaiwgy  w^Lii  wuLieo^uuuLU^  eiaatiu  ^iuuicuis  •  Lit 

the  procedure  consists  of  replacing  the  elastic  parameters  by 
expressions  involving  lineer  operators  that  contain  the  viscoelas¬ 
tic  properties  of  the  material  being  considered.  In  the  case  of 
steady-state  problems,  the  operator  expressions  reduce  to  alge¬ 
braic  functions  containing  the  frequency  as  well  as  viscoelastic 
properties  (see,  for  example.  Ref.  14).  For  the  sake  of  com¬ 
pleteness,  we  will  review  the  procedure,  following,  for  the  most 
part.  Ref.  14.  The  linear  strain  displacement  relations,  referred 
to  Cartesian  axes  (x^,  x^),  are 


*ij  ■  *  “j,i> 


i.j  -  1,  2,  3  . 


(3.4.1a) 


The  elastic  stress  strain  law  is 


oij  -  • 


(3.4.1b) 


The  equilibrium  equAt ion  is 


a2u. 


(3.4.2) 


In  Eqs.  (3.4.1),  and  in  what  follows,  i,j  ■  l,  2,  3;  repeated 
Indices  indicate  summation  over  all  the  possible  values  of  the 
index;  commas  denote  differentiation  with  respect  to  the  Car¬ 
tesian  coordinate  associated  with  the  index  following  the  comma, 
and  6^  is  the  Kronecker  delta  (B^j  »  0  when  i  j,  6^  ■  1 

when  i  «  j ) . 


The  stress  deviator, 
are  defined  as 


and  the  strain  deviator. 


* 


sij  "  °ij  "  ^6ij°kk  '  (3.4.3a) 

eij  "  eij  ’  *Bijekk  ‘  (3.4.3b) 
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The  creep  laws  are  written  as 


Ps-m  "  »  (3.4.4a) 


P',,idc-‘5’eklt'  »•*•«» 

I  I 

where  P,  Q,  P  ,  and  Q  are  linear  operators  with  constant 
coefficients  and  are  defined  as  follows: 
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(3.4.5a) 

(3.4.5b) 

(3.4.5c) 


(3.4.5d) 


Hie  constants  PQ,  P^,  , in  these  equations  are  determined 
on  the  basis  of  the  viscoelastic  model  selected. 


Since  the  present  analysis  utilizes  solutions  that  are  har¬ 
monic  in  time,  we  select  displacements  of  the  form 


(3.4.6a) 


where  barred  quantities  are  dependent  upon  space  coordinates 


only.  From  the  representation  (3. 
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6a)  f 

it  follows  that 

1 

L,j  +  ' 

(3.4.6b) 

'  eU 

“  ^6ijekk  ' 

(3. 4. be) 

'"atJ 

"  *5ij0kk  ' 

(3.4.6d) 

(3.4.6e) 

When  Eqs.  (3.4.6)  are  introduced  into  Eqs.  (3.4.4),  the 
creep  laws  become 


P(io))sij  -  Q(ia>)etj  , 


(3.4. /a) 


p'(ico)<j..  -  Q'(iai)e..  , 


(3.4.7b) 


where  now  the  symbols  P(ico),  Q(ico),  P^ico),  Q*(ico)  represent 
the  polynomials  in  (ico)  that  result  when  the  operations  defined 
by  (3.4.5)  are  carried  out,  that  is, 

F(ico)  “  Pq  +  ( ico) P ^  +  (ico)2p?  +  ...  +  (ico)pP^  ,  (3.4.8) 


If  we  define  the  following  rational  functions  of  (ico) , 


Q(ico) 

Yfl  -  (deviatoric  complex  modulus)  , 


Q*  (ico) 

Y  m  ■— (dilatational  complex  modulus)  , 
P  (ico) 


(3.4.9a) 


(3.4.9b) 


the  creep  laws  (3.4.7)  may  be  written  as 


sij  "  Vij  ' 


(3.4. 10a) 


°kk  "  Yv€kk  * 


(3.4. 10b) 


Then,  combining  Eqs.  (3.4.10)  with  the  definitions  of  the  deviator 
amplitudes  in  Eqs.  (3.4.6),  gives 


°ij  "  ^6ij*Yv  "  V€kk+  Ys€ij  * 


(3.4.11) 
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Comparison  of  Eq.  (3.4.11'  with  the  -"iastic  stress  strain 
Jaw,  Eq.  (3.4.1b),  Indicates  Chat  the  eia^'-ic  field  equations 
(and  corresponding  solutions)  may  be  emplo  /ed  Jn  o^oelestic 
problems  provided  that  the  elastic  modulii  X  and  p.  are  re¬ 
placed  by  the  corresponding  viscoelastic  modulii  X,  and  u  , 

V  v 

respectively,  where 


\  -  t(Yv  -  Ye)  ,  (3.4.12a) 

-  sYs  .  (3.4.12b) 


It  is  to  be  noted  that  X  and  p  are,  in  general,  complex, 
and  are  functions  of  the  frequency,  to,  and  the  parameters  that 
define  cha  viscoelastic  model. 

The  viscoelastic  equivalents  Xv  and  pv  must  replace  their 

elastic  counterparts  wherever  they  occur.  This  substitution  is 
necessary  not  only  where  X^  and  pm  appear  explicitly,  but  also 

in  those  parameters  in  which  they  are  implied,  namely,  or 

kp^  and  €.  tn  what  follows,  the  subscript  v  will  be  dropped, 
and  a  subscript  ir.  will  be  used  to  distinguish  the  medium  from 
the  liner. 

For  convenience  of  reference  we  rewrite  the  dependence  of 

these  parameters  on  X  paa  p  : 
r  m  m 


tsee  Eq.  (3.1.6c)],  which  may  be  put  into  the 

nondimens ional  form 


*>3  '  bko3 


[see  Eq.  (3.3.2e)  ]  , 


X  /(X  +  2p  )  [see  Eq.  (3.1.12c)], 

m  '  m  nr 


k  ,  ■  cu  /c  , 
p4  p  tra 


(3.4.13) 


or,  in  the  nondimens ional  form 


%  -  bV  ■ 
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We  now  j.nt::oduce  (3,4.12)  into  (3.4.13).  Recalling  that  Yv 
and  Yg  are  rational  functions  of  the  purs  imaginary  variable 
ioj,  and  that  rhp  r •equency  —  depends  u^uu  ehe  inae1*:  or  summa¬ 
tion  p,  (i.e.,  oj  •  pr/T)  we  wri^e 

Vtop>  -  Y5 p  *  Wvp5  ' 

(3.4. 14) 

W  "  Yrp>  +  1Ysp''  ' 

where 

Yd)  y(2)  YU)  v(2) 

vp  *  *vp  *  sp  ’  sp 

are  real,  and  the  subscript  p  is  used  to  indicate  dependence 
on  the  Index  of  summation.  In  addition,  we  let 

v(D  .  vvl)  +  2y(1) 

Yp  Yvp  ZYsp  * 

(3.4.15) 

Yn2)  "  Yv^  +  2Ys2)  ’ 
p  Vp  sp 

Then,  if  ws  introduce  (3.4.12)  into  (3.4.13),  and  take  account 
of  (3.4.14)  and  (3.4.15),  we  obtain  the  viscoelastic  equivalents 

of  k  ^  and  k  ,  in  nondimensional  form  (dropping  the  subscripts 
P3  P4 

p,  but  keeping  in  mind  the  dependence  on  p) ,  as  follows: 

q3  -  (bk3)  -  q,  -  iq3  ,  (3.4.16a) 

q4  -  (bk4)  -  q4  -  iq4  ,  (3.4.16b) 
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For  the  viscoelastic  counterparts  of  the  Lame  constant 
and  the  shear  modulus  we  obtain 


m 


X  +  iX  , 


(3 


__  'V 

u  -  p.  +  ip.  , 


.4 . 16c) 


,4.l6d) 


.4.  i6e) 


.4.16f) 


.4.17a) 


(3.4.17b) 


where 


A 


A 


M- 


M. 


i|Vl>  - 

Y^)  1 

l  vp 

sp  J 

*M2>  - 

V<2>] 

t  VP 

sp  J 

*£>  . 

sp 

}v<2>  . 

sp 

(3.4.17c) 


I  (3.4.17d) 


And  for  the  viscoelastic  counterpart  of  e. 


€  -  e  +  ie  , 

where 

—  A(A  +  2[x)  *h  A(A  +  2p.) 

®  _  —  9  ~  'v.o  « 

(A  +  2p-)Z  +  (A  +  2nr 

»v 

~  _  A(A  +  2p.)  ~  A(A  4-  2p.) 

(A  +  2p.)  +  (A  +  2|X)Z 


(3.4.18a) 


(3.4.18b) 


(3.4.18c) 


5.  VigggeiUaUc  Madel 

The  standard  linear  solid  viscoelastic  model  considered  in 
Ref.  15  and  shown  in  Fig.  4  was  selected  for  the  present  analysis. 
It  was  chosen  because  it  is  the  simplest  model  that  exhibits  the 
four  most  common  features  of  viscoelastic  behavior:  instantaneous 
elasticity,  creep,  stress  relaxation  and  creep  strain  recovery. 
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The  operator  representation  of  the  model  is 


_  'VmCdm 


T  i  + 


.LI 


m 


m 


ai  +  ^ 


/I  C  1  - \ 


m 


2  f 

7mctmLT2 _ 


Lt 


I 


(3.5.1b) 


where  the  notation  of  Ref.  15  has  been  changed,  and  the  symbols 
employed  above  have  the  following  meaning: 


velocity  of  propagation  of  high  frequency 
4  dilatational  waves 

c  velocity  of  propagation  of  high  frequency 

distortional  waves 

t  strain  recovery  tine  In  uniaxial  strain  case 


.L 

T2 


recovery  time  for  shear  strain 


relaxation  time  for  normal  stress  in  uniaxial 
strain  case 


relaxation  time  for  shear  stress. 


The  choice  of  the  specific  model  defined  by  (3.5.1)  now 
enables  us  to  evaluate  the  operators  represented  in  general  form 
in  Eq,  (3.4.5)  and  associated  with  the  stress  strain  relations 
in  (3.4.4; . 
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From  (3.4.12b)  and  (3.4.9a), 
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The  present  analysis  is  concerned  with  steady-state  solutions 
i%t 

in  which  the  factor  e  v  is  present.  Therefore,  when  the 
operators  Q,  P,  Q1 ,  and  P1,  as  defined  by  (3.5.2)  and  (3.5.3), 
iAj.e  uacu  in  uie  li cep  lawa  (3*4.4),  they  generate  the  polynomials 
Q  (ioDp) ,  P(icOp),  Q'(iu>p),  p'(io)p)  indicated  in  (3.4.8).  These 

polynomials  now  take  the  form: 


Q 


2y  c 2 

‘m  tm 


+  loo 

P 


> 


P  -  +  ia)p  * 

o'  -  %  +  -  Vp 

it  i  2 

P  -  Pn  +  io>  P.  -  03  , 

0  pi  p  * 


where 


(3.5.4) 


% 


dmTlfi2 


i 


0, L 


'm 


3<4(t1  +  n2>  -  4ctm(T2  +  ni> 


I 


0  "  nin2  ' 

Pi  ■  J?i  +  • 


The  deviatoric  and  dilatational  complex  moduli!  Ys  and  Yv» 
Eq.  (3.4.9),  may  be  written  in  terms  of  the  basic  properties  of 
the  viscoelastic  model  by  utilizing  the  polynomials  obtained 
above. 
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i 


By  using  (3.5.4),  the  devlatoric  modulus,  Y 

8 


Q(to^) 

Pa<V 


may  be  written  in  terms  of  real  and  imaginary  parts  as 


Y 


s 


+  iY 


(2) 

'80 


9 


where 


,( D 

sp 


2y  c  l 
'm  tm 


T2  (  .  2  /cdtt/n2b\' 

S*(pr) 


1  +  <p*>2 


(3.5.5) 


„C2>  . 

sp 


2VtmpT 


i  -  la 
1 

c,  /to-b 

U.  T/b  ) 
-  dm 

1  +  w*  (^i^) 


Q  C^^p) 

Similarly,  the  dilate tional  complex  modulus,  Y„  ■  , ,  , 

v  P  (i<0 
P 

may  be  written  in  terms  of  its  real  and  imaginary  parte  as 

Y  -  Y(1)  +  iYC2) 

Yv  Yvp  +  1Yvp  * 

where 

V<U  .  vUK.'  Q2<p^2][p0  ~  (P^)2]  +  Wpir)2 


'vp 


Pq  "  (P’O2  +  P2(ptt)2 


,(2)  .  ^cL(PT>{[P0  ~  <Pt)2]%  ~  K  '  Vplr)2 


il 


(3.5.6) 


vp 


pq  *  (p^)2  +  p2  (p,r) 2 
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where,  froraEq.  (3.5.4), 


i  9  or 

c,  T/b  .  ~  .ft. 


A  *1 


a  .  Sli-  .  (?*? )  A\|s  II  .  «/2t» 

%  2  VcJ_/a,w  InJra.  4U .  ; 

y  c ,  uui  a  z  i  i  am 

m  am  l 


h 

n 

2 


\  ' 


V 

7 

m  dm 


J3 

■Wn2bj  I 


(|)  (£)  ^  b 


\C  ,  / 

dm 


Z2  +  £l 

Q2  a2 


}■ 


y  Cj 

m  dm 


Qo  -  ~ V*  -  3  - 


<£)  ’ 
dm 


(3.5.7) 


;  _  p't2  .  (°*^'bA  A) 

F0  P°T  Cdm^2b^  -2  ’ 


_  ,  c ,  l/b  ft.^ 

*i-V-(;£*SE>  (*♦«$>  • 


Note  that,  by  the  way  Eqs.  (3.5.5),  (3.5.6),  and  (3.5.7)  are 
written,  we  have  introduced  in  the  formulation  of  the  problem 
the  following  non dimensional  ratios  to  represent  the  visco¬ 
elastic  properties  of  the  medium: 


shear  relaxation  time  to  uniaxial  strain 
relaxation  time 


ll 

ft. 


stress  relaxation  timi  to  strain  recovery  time 
for  uniaxial  strain 


h 

fto 


shear  stress  relaxation  time  to  shear  strain 
relaxation  time 


cdm 

r‘2b 


shear  relaxation  time  to  half  transit  time. 
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The  ratios  and  als°  represent  the  ratios  of  the 

relaxed  to  unrelaxed  elastic  properties  of  the  medium,  namely, 
in  the  notation  introduced  in  Eq.  (3.4.17), 


T, 

.4 


(2u  + 


(2|i 


J3. 


m 


A  ^ 
nr  co 

+  A  )  P 

m'co 


-»  0 
-»  w 


(n 


CO 


-»  0 
« 


(3.5.8) 


This  can  be  verified  by  letting  p  ->  0  and  «  in  Eqs.  (3.5.5) 
and  (3.5.6)  and  referring  back  to  Eq.  (3.4.17).  Note  that  for 
finite  a^,  the  ratios  in  (3.5.8)  can  also  be  obtained  by  letting 

c  ,  /fl„b,  instead  of  cu  ,  approach  0  and  oo. 
dm  i.  p 


6 .  Stresses  and  Displacements  Applied  to  the  Cylinder  by  an 
Incident  Dilatational  Plane  Wave  in..a  VlscQ&Lastlc  Mfidium 

In  this  section,  we  consider  the  elastic  cylinder  to  be  enr 
bedded  in  a  viscoelastic  medium,  and  obtain  the  stresses  and  dis¬ 
placements  applied  to  the  cylinder,  by  an  incident  dilatational 
plane  wave.  In  doing  this  we  use  the  form  of  the  equations  de¬ 
veloped  in  Section  III.l  for  the  case  of  an  elastic  embedding 
medium,  and  invoke  the  correspondence  principle  discussed  in 
Section  III. 4.  That  is,  we  take  the  results  obtained  for  an 
elastic  medium  and  replace  the  parameters  involving  elastic  modulii, 
by  their  viscoelastic  equivalents. 

For  convenience  of  reference,  we  rewrite  the  pertinent  results 
obtained  in  Section  III.l.  From  Eqs.  (3.1.13)  the  normal  and 
shearing  stresses  applied  to  the  cylinder  boundary  by  the  incident 
plane  dilatational  wave  are  given  in  polar  coordinates  by 

oo 

_  iCDt 

V  0  e  p  ,  (3.6. la) 

^  r,p 

P-1 

_  “  io>  t 

<£»  -  -  l  v,p  *  p  -  <3-6'lb) 

P-1 
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in  which 


2  2  iqo(cos  0-1) 

o..  «  (-ia  ) (cos^e  +  e  sin  0)  e 

*•  »F  F 


M.6  .  Ir'l 


rS 


,P  ■  {‘UP>  (^T*)  sln  29  e 


iq^Ccos  e-l) 


(3. 6. Id) 


where  the  viscoelastic  parameters  and  e  are  given  by 
Eqs.  (3. A. 16)  and  (3. A. 18)  in  the  form 


q3  ‘  q3  "  iq3  ' 


(3.6.2a) 


e  -  e  +  ie 


(3.6.2b) 


As  in  previous  cases*  we  expand  the  boundary  tractions  in 
the  Courier  series 


00  00 


id)  t 

4  -  l  I  aPV c08  "9  e  p 

p«l  n»0 


(3.6.3a) 


Y  y  a  B  sin  n 9 

L  L  p  pn 


i(i>  t 


r0  L  L  p  pn 

P"1  n«l 


e 


(3.6.3b) 


and,  comparing  (3.6.3)  with  (3.6.1),  we  obtain  the  coefficients 

A  .  B  as 
pn  pn 


& 


A  -  -  U 

pn  -rr 


2  2  iq3(cos  e-l) 

(cos  6  +  e  sin  e)e  cos  ne  de 

n  >  1 


(3.6.Aa) 
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Adq  -  -  “  (cos20  +  e  si.,2,-)  e 


2  iq^ (cos  h-L) 
i')  e  d-7  , 


(3.6.4b) 


i (1  -  €)  iq3(cos  b-1) 

B  „  ■  - —  —  e  sin  20  sin  n£  d 9  . 

pn  ir 


(3.6.4c) 


These  expressions  are  evaluated  in  terms  of  Bessel  functions 
of  a  complex  argument  in  Appendix  B.  As  given  in  (3.6.4),  they 
are  identical  in  form  to  their  elastic  counterparts,  (3.1.15), 
but  of  course  are  different  in  value. 

A  minor  variation  is  necessary  in  the  development  of  the 
displacement  representation.  From  Eq.  (3.1.17), 


ice  t 

l  Vp  ‘  P  ' 


(3.6.5a) 


_  lew  t 

-  V  „  *  P 


“e  ~  I  ue,p  -  ' 

p*l 

in  which  we  take  (3.1.17e  and  f)  as  the  form  appropriate  for 
the  viscoelastic  representation  of  ur  u0  ,  Thus, 


(3.6.5b) 


aP  t  cdir^3\  iq3(ccs  9-1) 

u  - - - -  ■'a\*)cos  0  e  i 

r,p  y  c  .  0)  \  Oi  b  / 

dm  p  p 


(3.6.5c) 


±P  ,  _i<l3(cos  '’'l> 


-x—  ( -a"“)sin  0  e 

c  ,  o)  \  cu  b  ) 


(3.6. 5d) 


m  dm  p  p 


where  we  note  again  that  q3  is  a  viscoelastic  parameter,  de¬ 
fined  by  (3.6.2a)  and  the  symbol  cjm  is  the  speed  of  propagation 
of  dilatational  waves  in  an  elastic  medium  that  has  the  unrelaxed 
elastic  properties  of  the  viscoelastic  material. 
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R  R 

Expanding  ch<j  displacement u  ind  a 


n  me  t  c.it  ief 


ser iet 


00  CO 

B  v  f 

u  .  -  -  >  <  ) 


)  <  )  -  C  cos  n0l  e 

L  |  L  y  c.  at  pn  | 

v.  m  d:  i  n  r  J 


n  la)  t 


p«l  n-0 


r  w  j 

ro  du  p 


(3.6.6a) 


B 


co  Oo 

Fly—  - 

L  1  L  v  c,  a)  pn 
.  v-  i  m  dm  p  v 
p-1  n-1  v 


1  ico  t 

D  s in  n0 t  e  ^ 


(3.6.6b) 


and  comparing  (3.6.6)  with  (3.6.5),  we  obtain 

& 


0 

pn  tt  \  co  b  / 

P 


0 


iq3(cos  0-1) 

e  cos  0  cos  nd  d0  , 

n  >  1 


C  .  1  fcdmq3\ 

up0  TT  \  tD  b  / 

P 


iq3(cos  0-1) 

e  cos  0  d0  , 


(3.6.7a) 


(3.6.7b) 


Dpn  " 


1  /fdafo 

7T  \  £U  b  / 

P 


“o 


iq3(cos  0-1) 


sin  0  sin  n0  d0  . 


(3.6.7c) 


The  form  of  the  quantities  in  (3.6,7)  differs  from  the 
counterparts  in  Eq.  (3.1.19)  only  by  the  presence  of  the  factor 
(Cdnj^/^pb)  •  When  the  parameter  q3  applies  to  an  elastic 

medium,  this  factor  is  unity.  The  explicit  representation  of 
(3.6.7)  in  terms  of  Bessel  functions  of  a  complex  argument  is 
given  in  Appendix  B. 
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(3.6.4b) 


P« 


,  2..  .  2  .  ^3U°S  "1}. 

(cos  +  f  sin  )  e  (1  , 


-7r 


i(l  -  «) 


pn  7T 


iq^ (cos  o-l) 

e  sin  2 <j  sin  nu  d 0 


(3.6.4c) 


These  expressions  are  evaluated  in  terms  of  Bessel  functions 
of  a  complex  argument  in  Appendix  B.  As  given  in  (3.6.4),  they 
are  identical  in  form  to  their  elastic  counterparts,  (3.1.15), 
but  of  course  are  different  in  value. 


A  minor  variation  is  necessary  in  the  development  of  the 
displacement  representation.  From  Eq.  (3.1.17), 


u 


r,P 


io>  t 

e  F  , 


(3.6.5a) 


n  ico  t 

u0  "  I  u0,p  e  P  '  (3.6.5b) 

P-1 

in  which  we  take  (3.1.17e  and  f)  as  the  form  appropriate  for 
the  viscoelastic  representation  of  ur  ,  nQ  .  Thus, 


aT 


y  c.  to  V  ou  b  ) 

•m  dm  p  p 


'O. 


iq.j(ccs  0-1) 
e  > 


<3.6. 3c) 


y  c  ,  o> 

rm  diu  p 


\  u>  b  l 

P 


sin  0 


iq^cos  0-1) 
e 


(3.6 . 5d) 


where  we  note  again  that  q3  is  a  viscoelastic  parameter,  de¬ 
fined  by  (3.6.2a)  and  the  symbol  c^  is  the  speed  of  propagation 
of  dilatational  waves  in  an  elastic  medium  that  has  the  unrelaxed 
elastic  properties  of  the  viscoelastic  material. 
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7 •  Solution  to  the  Wave  Equation  In  d  Vlscuelast 1c  Medium 


In  Section  III. 3,  the  solution  to  the  wave  equation  In  a 
surroundine  elastic  medium  w »a  rnn af-rnri-nrf  displaccnt  nt 

potentials.  Ihe  displacements  and  stresses  associated  with 
this  solution  are  given  by  Eqs.  (3.3.2)  and  (3.3.4).  Bv  utiliz¬ 
ing  the  correspondence  principle  (Section  III. 4),  we  nay  obtain 
the  solution  in  a  viscoelastic  medium  from  the  solution  in  the 
elastic  medium.  The  modulii  appearing  in  the  elastic  solution 
are  simply  replaced  by  their  viscoelastic  equivalents.  Thus, 
the  desired  solution  in  the  viscoelastic  medium,  that  is, 
stresses  and  displacements,  are  given  by  Eqs.  (3.3.2)  and 
(3.3.4),  provided  that  the  parameters  appearing  in  these  equa¬ 
tions  are  defined  as  follows: 


k  -  ■  k  - 
p3  p3 

'V 

-  ik 

1 

/V 

A  -A 
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m 

m  V 

ta 

+  ip. 

(3.7.1) 


For  convenience  of  reference,  the  stresses  and  displacements 
given  by  Eqs.  (3.3.2)  and  (3.3.4)  are  rewritten  here,  with 
Eq.  (3.7.1)  now  applying. 
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(3.7.2b) 
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(3. 7. 2d) 


(3.7.3a) 


(3.7.3b) 
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It  1b  observed  that  the  viscoelastic  parameters  and  ^ 

necessitate  the  computation  of  Bessel  functions  of  a  complex 
argument. 


8.  Superposition  of  Solutions  to  Enforce  the  Boundary  Conditions 

In  the  preceding  sections,  we  have  presented  contributions 
to  the  complete  solution.  These  consisted  of;  1)  stresses  and 
displacements  associated  with  the  incident  wave;  2)  the  solution 
to  the  wave  equation  in  the  elastic  cylinder;  and  3)  the  solution 
to  the  wave  equation  in  the  surrounding  medium.  These  contribu¬ 
tions  must  be  combined  to  enforce  the  required  conditions  at  the 
cylinder  boundaries,  namely,  that;  1)  the  inner  boundary  be  trac¬ 
tion  free;  and  2)  at  the  outer  boundary  the  cylinder  displace¬ 
ments  and  appropriate  stress  components  be  continuous  with  those 
of  the  surrounding  medium. 


Recalling  the  introduction  of  the  nondimen sional  radial 
coordinate  p  ■  r/b  [Eq.  (3.2.9c)  ],  and  letting  p  ■  a/b,  we 
obtain  the  boundary  conditions  at  p  ■  £  (inner  boundary  of  the 
cylinder)  in  the  form 
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o°  (;•)  +  o's  (_)  -  u  , 
rr  '  rrx  '  ’ 


d  ^  ^  8 

°r*  +  vro^  "  C  • 


,(3. 8. la) 


At  p  «  L  (outer  boundary  of  the  cylinder). 


<&.U)  +  o’  (l)  -  of“(l)  -  o“(l)  - 


rr 


rr 


rr 


rr 


rr 


o%a)  +  cj5(l)  -  $(1)  -  0*“(1)  -  4  , 


u£(l)  +  uj-(l)  -  uj°(l)  -  u£m(l) 


-uB 


(3.8.1b) 


u?(D  +  u!(i)  -  u^a)  -  u!mci)  -  uB 


where:  1)  the  superscripts  d  and  s  are  affixed  to  quantities 
s terming  from  the  dilatational  and  shear  solutions,  respectively, 
to  the  wave  equation  in  the  cylinder  (see  Section  III. 2);  2)  the 

superscripts  dm  and  sm  are  appended  to  stresses  and  displace¬ 
ments  resulting  from  the  dilatational  and  shear  solutions,  respec¬ 
tively,  to  the  wave  equation  in  the  medium  (for  an  elastic  medium, 
these  quantities  are  given  in  Section  III. 3;  for  a  viscoelastic 
medium,  they  are  given  in  Section  III. 7);  and  3)  the  superscripts 
B  identify  the  stresses  and  displacements  associated  with  the 
incident  wave  at  the  outer  boundary  of  the  cylinder  (for  a  wave 
traveling  through  an  elastic  medium,  these  quantities  are  given 
in  Section  III.l;  when  the  wave  propagates  through  a  viscoelastic 
medium,  they  are  given  in  Section  III. 6) . 


In  the  following  development,  advantage  will  be  taker*  of  the 
fact  that  by  an  appropriate  choice  of  the  values  of  the  visco¬ 
elastic  parameters,  the  viscoelastic  solutions  in  the  medium  re¬ 
duce  to  elastic  solutions.  Therefore,  quantities  labeled  with 
the  superscripts  B,  dm,  and  sm  will  be  written  as  they  are 
given  in  Sections  III. 6  and  III. 7.  Thus,  the  solution  presented 
will  apply  to  an  elastic  cylinder  embedded  in  a  viscoelastic 
medium,  and  as  a  limiting  case  wiil  contain  the  solution  for  an 
elastic  cylinder  embedded  in  an  elastic  medium. 


42 


Ai'ter  the  terms  entering  into  Eq.  (3.6.1)  are  computed  from 
the  appropriate  expressions,  and  the  common  trigonometric  and 
exponential  factors  are  divided  out,  the  resulting  nonhomogeneous 
system  of  linear  equations  in  the  unknown  coefficients  of  super- 

a  4  4  *7  mo  t.rv  <  f-  A  .#*»  mn  I  v  n  r. 

r  ~j  •  Pn  "  . . . . . . 


[CtJ/pn 


jz  \  -  |d  } 

1  \  i»pnj 

p  a  1,  2,  3,  • . . 


i,j  «  1,  2,  6 

1»  2,  ...  , 


(3.8.2) 


where  C,, _ stems  from  the  coefficient  of  Z,  in  the  ltl1 

i J / Pn  j i pn 

equation  of  Eqs.  (3.8.1);  D.  results  from  the  term  on  the 

th  LiPn 

right-hand  side  of  the  i  equation  of  Eqs.  (3.8.1);  and 
(Z^pnJ/  the  solution  vector,  is  a  6x1  column  matrix,  whose 

complex  elements  are  the  unknown  coefficients  of  superposition. 
It  should  he  noted  that  there  is  a  different  matrix  equation 
(3.8.2)  for  each  pair  of  values  of  p  and  n. 


The  element#  C, .  __  and  D,  are  shown  in  detail  in 
ijjPn  i/pn 

Appendices  A  and  3. 


Ihe  system  (3.8.2)  may  be  solved  by  m  trix  inversion  for  the 

solution  vector  (Z.  ‘  ).  Having  solved  for  [Z.  )  for  a  suffi- 

jipn  J,pn 

cient  number  of  values  of  p  and  n  to  ensure  adequate  conver¬ 
gence,  we  may  calculate  to  the  desired  accuracy  the  stresst?  and 
displacements  in  the  cylinder  as; 


rr 


aa  +  os 

rr  rr 


d  .  s 
rtf  rtf  c  rtf 


'99 


_d  .  £ 

aee  +  aGG 


, (3. 8,3a) 


ur  "  ur  +  ur  '  u9  "  utf  +  u*  »  (3.8.3b) 

where  the  terms  appearing  on  the  right-hand  side  of  the  above 
expressions  are  given  in  Eqs.  (3.2.10)  and  (3.2.17). 
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SECTION  IV 


SOLUTION  FOR  AN  IMClPEMI^ISlQmOMAL  PLANE  WAVE 


IH  €5  /^11  a  a  4  r>r\  rsG  A  4  V  f  ^ 

- r--  Tf - T.4>  . 


This  section  is  concerned  with  the  cylinder  response  to  an 
incident  shear  wave.  As  discusseo  in  Section  III. 4,  the  visco¬ 
elastic  analysis  may  be  obtained  from  the  corresponding  elastic 
analysis  by  using  the  correspondence  principle,  that  is,  by  re¬ 
placing  the  elastic  modulil  by  appropriate  viscoelastic  moduli!. 
Conversely,  the  elastic  case  can  be  obtained  as  a  special  case 
of  the  viscoelastic  formulation  by  a  proper  choice  of  the  value 
of  certain  parameters.  Therefore,  in  treating  the  response  to  an 
incoming  shear  wave,  we  shall  regard  the  cylinder  as  first  em¬ 
bedded  in  an  elastic  medium.  The  viscoexastic  analysis  will  then 
be  obtained  by  employing  the  correspondence  principle.  The  numer* 
cal  results  will  be  computed  from  the  viscoelastic  formulation  for 
both  the  elastic  and  y is coelastic  case,  treating  the  elastic  prob¬ 
lem  as  a  limiting  case  of  the  viscoelastic  problem. 

The  first  section  In  this  chapter  .deals  with  the  .contribution 
of  the  free  field  quantities,  that  is,  the  stresses  and  displace¬ 
ments  applied  to  the  outer  boundary  of  the  cylinder  by  the  inci¬ 
dent  shear  wave.  This  contribution  will  be  derived  independently 
of  any  of  the  results  presented  in  Section  III.  ttie  remaining 
contributions  to  the  solution,  namely,  the  solution  to  the  wave 
equation  in  the  cylinder  and  the  scattered  response  in  the  sur¬ 
rounding  medium,  will  be  obtained  by  making  minor  changes  in  the 
corresponding  solution  for  an  incoming  dllatational  wave.  As 
outlined  in  Section  III, 2,  these  alterations  involve  interchanging 
«i-i  tS*  and  cos  n 0,  and  setting  K  •»  -  n.  It  will  be  recalled 
that,  the  TT  are  those  integers  stemming  from  differentiation  with 
respect  to  0 . 


2.  Stresses  and  Displacements  Applied  to  the  Cylinder  by  an 
Incident  Distortlonal  Plane  Wave 

To  avoid  confusion  with  the  response  to  ar  incoming  dilata- 
tional  wave,  we  shall  identify  corresponding  quantities  associated 
with  the  response  to  an  incoming  shear  wave  by  using  notation 
that  differs  slightly  from  that  employed  in  Section  III. 
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The  Cartesian  coordinate  system  in  which  the  incident  plane 
cistorticnal  wavj  is  described,  ae  we.1 1  as  some  of  the  nota-von 
used,  is  shown  in  Fig.  3b.  The  disturbance  is  assumed  to  be 
traveling  in  ; he  negative  x-directicn.  Ihe  propagation  of  the 
v.’avc  depends  only  on  tlic  apaoc  ouuiuinaue  x,  ana  the  time  coordi¬ 
nate  t. 

We  may  obtain  the  governing  equations  of  mtion  from  a  scalar 
potential,  T,  Ref.  9,  where  T  -  ¥(£,  t)  and  is  a  solution  of 


ix  -  7  f  •  0 

m  ox2  ra  dt2 


(4.2.1) 


The  nonvanishing  Cartesian  displacement  and  t.  ress  v'anti" 

ties  v_  and  t_-  are  given  by 
v  xy  "  3 


v— 

y 


.  $1 
dx  ’ 


(4.2.2a) 


T--  *  -  u  -  (4.2.2b) 

xy  ®  s-2 

As  developed  in  Ref,  9,  Eqs.  (4.2.1)  and  (4.2.2)  apply  to 
an  elastic  medium,  and  pm  is  the  elastic  shear  modulus.  How¬ 
ever,,  on  thu  basis  of  the  correspondence  principle,  these  equa¬ 
tions  may  also  be  applied  to  a  viscoelastic  medium  if  is 

regarded  as  the  viscoelastic  parameter  defined  in  Section  u,4 
by  Eqs.  (3.4,17b)  and  (3.4.17d'1  and  given  in  terms  of  the  visco¬ 
elastic  model  representation  used  in  the  present  analysis  by 
Eq.  (3.5.5). 

Steady  state  solutions  of  (4.2.1)  may  be  taken  in  the  form 


T  -  T  (5c)  e 


ico  t 


(4.2.3a) 


where  the  frequency  a>  is  given  by 


oo 


p7T 

T  ' 


(4.2.3b) 
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wi  Lh  p  a  posit iv-  im.fpcr,  and  X  an  arbitrary  time  interval. 
Thus,  the  wave  equation  (4.2.1)  becomes  the  HeLmhoLz  equation 


where 


,2... 

~  jP 

a*2 


k2  T 
p4  p 


-  0 


2 

tip 


(4.2.4a) 


(4.2.4b) 


In  the  case  of  a  viscoelastic  medium, 

parameter  defined  in  Ea.  (3.4.16).  Fcr  an 

reduces  to 


is  the  complex 


elastic  medium. 


(4.2.4c) 


where  .c  is  the  velocity  of  propagation  of  -die  tort  ional 
waves  in  that  medium,  with 


(4.2. 4d) 


Steady-state  solutions  of  (4-2,1)  may  be  constructed  from 
the  solutions  to  (4.2.4a),  by  sing  (4.2.3a).  Thus  T(x,t) 
may  be  taken  as 


T(x,t)  -  £  A  e 


i[kp4(x-b)  +  wpt  J 


(4.2.5) 


P-1 


We  now  represent  the  incoming  shear  wave  in  the  form 


x__  ■  )  a  sin  co 

xy  4.  p  P 

p-1 


co 


(4.2.6) 
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and  proceed  to  relate  the  coefficients,  ap,  in  this  equation  to 
the  coefficients.  A_.  in  Ea.  (4.2.  SI.  Tn  th<»  olooHr  ^»?e 

r  ■ 

ojp/kp^  «  ctm,  t^ie  un^orra  speed  of  propagation  of  distortionai 

waves,  in  the  viscoelastic  case  the  parameter  co^/kp^  represents 

the  limiting  value  of  the  velocity  of  propagation  when  the  fre¬ 
quency  tends  to  infinity.  At  x  -  b, 

t  *  Y  a  sin  to  t  .  (4.2.7) 

Sy  L*  P  P 

p-1 


The  Fourier  coefficients  3p  define  the  wave  form  over  the 

arbitrary  half  period,  T.  In  the  viscoelastic  case,  since 
attenuation  occurs,  the  wave  form  is  altered  as  it  progresses 
through  the  medium. 

An  expression  for  t_.  identical  with  (4.2.6)  is  given  by 

xy 

the  real  part  of 


T„. 

xy 


“  i[k  4(x-b)  +  co  t] 

l  <-%> B 

p-1 


(4.2  ) 


From  (4.2.2b)  and  (4.2.5), 
00 

2 


t  -  Y  |i  k  ,  A  e 
5y  L  rn  p4  p 


iCkp4(x-b)  +  o)pt] 


p-1 


(4.2.9) 


Comparing  (4.2.8)  and  (4.2.9),  we  obtain 

ia„ 


^m^p4 


(4.2.10) 


From  (4.2.2a)  and  (4.2.5),  and  by  noting  (4.2.10),  we  get 


v 

y 


oo 


i. 

p-i 


aP 

p  k  , 
ra  p4 


i[k  4(x-b) 
6  V 


+  CUpt] 
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It  will  be  convenient  for  the  computational  effort  to  rewrite 

2 

the  above  expression.  By  (4.2.4b),  p  k  ,  -  (y  ou  /k  .)  and  from 

m  p4  % 'm  p  p4 

the  definitions  of  co  [Eq.  (4.2.3b)]  and  q,  [Eq.  (3. 4. lb'll. 

?  •+ 

we  may  rewrite  p  k  .  as 
m  p4 


1 _  „  _ l _  /  Ctmq4\ 

k  ,  v  c,  a)  \  ra  b  /  ‘ 


Thus, 


I 


m  p4  m  ttn  p 


an  q.  i[k  ,  (x-b)  +  to  t  ] 

_ E _  e  p4  P 

y  c  co  \  co  b  / 

^  m  tm  p  p 


(4.2.11) 


In  polar  coordinates  the  stresses  and  displacements  associa¬ 
ted  with  the  incoming  wave  are  given  by 

T®  -  T_.  sin  26  ,  -  T  cos  26  , 

rr  xy  rd 


v„  sin  6  ,  v*  -  V_  COS  9  . 

V  &  V 


Substituting  (4.2.8)  into  (4.2.12),  and  setting  x  ■  b  cos  0, 
we  obtain  the  following  expressions  for  the  stresses  applied  to 
the  boundary  by  the  incident  wave: 


TrS  "  L  Tr£, 
P-1 
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where 


_  iq, (cos  0-1) 

Tr  n  -  (“iajflin  2 9  e  , 

»  »  * 

_  iq,(cos  0-1) 

Tr?,P  "  (*iap)(c0S  29)  6 


(i.2.  lit.' 


(4.2. 14d) 


with 


qA  -  bkp4  * 


(4.2. 14e) 


is  defined  in  Eq.  (3.4.16)  for  the  viscoelastic  case.  In 
the  elastic  case  it  simplifies  to 


baip 

:.tm 


q4  "  c  * 


The  Fourier  expansions  of  Eq.  (4.2.14)  in  the  circumferen¬ 
tial  coordinate  £  are 


Trr  *  I  l  *pApn  sln  n?  e  P  * 


ico  t 


rr  L  Li  p  pn 
p«l  n«l 


(4.2.15a) 


_  ico  t 

,  a  B  cos  n©  e  " 
L  p  pn 

p»l  n-0 


I  l 


(4.2.15b) 


The  complex  Fourier  coefficients  Apn  and  Bpn  are  obtained 
by  comparing  (4.2.15)  with  (4.2, 14) >  and  are 


& 


-  - 


pn  t r 


iq, (cos  0-1)  _  _  _ 

e  sin  20  sin  n0  d©  , 


(4.2.16a) 
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n  >  1  (4.2.16b) 


i 


pn 


u 

7 r 


iq.  (cos  9-1)  0T  “ 

n4 '  cos  2-  cos  n9  d'J  , 


JT 


po 


i 

t r 


iq^ (cos  9-i)  _  _ 

e  cos  29  d9  . 


(4.2.16c) 


The  representation  of  these  coefficients  in  terms  of  Bessel 
functions  of  a  complex  argument  is  given  in  Appendix  B. 

If  we  introduce  (4.2.11)  into  (4.2.13)  and  set  x  ■  b  cos  9, 
we  obtain  the  following  representation  for  the  displacements 

B  .  B 
vr  and  vg : 


B 


l 

P-1 


icu  t 

v,p  •  '  * 


(4.2.17a) 


B 

ve 


- 1 

p-i 


i»  t 

ve,P  6 


(4.2.17b) 


where 


a. 


r»P 


(sf)  % 


Vt 


7  ct  0) 

'm  tm  p  p 


e 


iq4(cos  ?-I) 


sin  9  , 


0,p  -  K  ^4  cos  9  e 

yv  rm  tm  p  p 


_  iq.(cos  e-i) 

a  *  H 


(4.2.17c) 


(4.2. 17d) 


To  facilitate  the  application  of  conditions  of  continuity 

of  displacement  at  the  outer  boundary  of  the  cylinder,  we  expand 

B  B  ** 

v  and  V?  in  Fourier  series  in  9.  Thus, 
r  9  ' 
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00 


oo 


B 


-  )  V  7rr—r 

",  'nTtm  p 
p-1  n-1  r 


_  ito  t 

C _ si.i  nfl  e  ** 

p.« 


00  oo 


B 

v* 


-  I  l 


io>  t 


0  ^  L  7  c.  0)  pn 

.Jn  'i  tin  p  v 
p*l  n»0  r 


D  cos  n0  e 


(4.2.18b) 


Cougar ing  (4.2.18)  with  (4.2.17),  we  obtain  for  the  coeffi- 

cients  C  and  D  . 

pn  pn 


.1 T 


^pn  "  i r  \o>  b/  q4 
P 


iq.(cos  0“1)  «  _ 

e  sin  0  gin  nd  d0  , 


(4.2.19a) 


V  ‘  ‘  q4 

P 


iir 


iq,(cos  0-1) 

a  cos  0  cos  n0  d0  , 

n  >  1 


(4.2.19b) 


«7T 


]}  m  ^  q 

pO  ir  \tu  b/ 


iq. (cos  0-1)  ■_  _ 

e  *  cos  0  d0  . 


(4.2.19c) 


These  coefficients  are  represented  in  terms  of  Bessel  func¬ 
tions  of  a  complex  argument  in  Appendix  B. 


3.  Solution  to  the  Wave  Equation  in  the  Elastic  Cylinder 

We  have  already  discussed  the  manner  in  which  it  is  possible 
to  apply  the  analysis  derived  for  the  response  to  an  incident 
dllatational  wave  to  the  corresponding  analysis  for  the  response 
to  an  incident  shear  wave  (see  Section  IV. 1  and  the  comments  pre¬ 
ceding  Eq.  (3.2.8)  in  Section  III. 2).  Thus,  if  in  the  results 
obtained  in  Section  III. 2,  we  Interchange  sin  n0  and  cos  nd 
and  regard  !T  as  equal  to  -n,  we  may  write  the  contribution, 
associated  with  waves  of  dilatation  in  the  cylinder,  in  response 
to  an  incident  shear  wave,  as  follows: 
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d  1 

i  «  — 


b  L  L  qijzi,i 

n-1  n-1 


pnLq.p  ^n^ql^  ^n+l^ql^ 


+  Z 


2, pi 


tj.  W>  '  VlM  ]} 


sin  n@  e  ^ 


(4.3.1a) 


+  £  Apl  sin  I  , 
P-1 


VS  ‘  t  i  L  P  KjaW)  +  Z2,PnW>}“8  n? 
p-1  n-0 


+  I  v  co* 8 ' 


id)  t 

>  P 


(4.3.1b) 


P-1 


00  00 
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p-1  n-1 


-l#pn 


l)K(v> 
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(4.3.1c) 


{x  ■  •  i)K(v> 
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_  la)  t 
sin  n©  e  ^  , 


03  ^ 

-"‘hi  1  ql  Zl,pr.  {A  + 

UU  0-1 


p»i  n«l 


2tJ.n(n  *  1) 

(q,p)2 
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(qLp)2 
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P)  (4.3. Id) 


ia>  c 
.  P 


+  Yn+l(qlp)  |  flln  nS  e 


l  l  J2{Zl,im[V'WV)  •  (n  ' 


r.Q  b  %  p 

p»l  n»0 


r  :  »»  .  io>  t 

+  Z2,pnj  <V>  WV>*  <"  ‘  ]““>•  n*  «  ° 


(4.3. le) 


where  the  superscript  d  has  been  used  to  designate  the  contribu¬ 
tion  from  the  dilatations!  wave  solution  in  the  cylinder. 

The  contribution  from  the  shear  wave  solution  in  the  cylinder 
is  obtained  from  Section  III. 2,  in  a  similar  manner,  and  is  given 
as  follows:  ’  ...  : 

'  „  '  ■  ■  . ■  .  :•'•  ••  ■ 

oo  oo  .  •  ••;vr 

vr  ’  -b  I  I  ?{Z3,pnW> 
p«l  n-1 

1  (4.3.2a) 

v'  ■■■ '  .  _  tot  4  _ 

+  Z4,pn(q2p)Yn(q2p) }sin  e  P  *  I  Bpl  8tn  ®  * 
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b  l  I  ^KpJ  .^7  W> 

P“i  n— C 
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(4.3.2b) 
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4.  Scattered  Wave  Solution  In  the  Surrounding  Medium 


In  this  section,  we  write  the  stresses  and  displacements 
associated  with  the  scattered  wave  solution  in  the  surrounding 
moHlum  And  which  contribute  to  the  resDonse  of  the  cylinder  to 
an  incoming  shear  wave.  In  doing  this  we  use  the  results  of 
Section  III. 3,  with  sin  n 6  and  cos  n 6  interchanged  and 

n  -  -  n.  Thus, 
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In  the  preceding  expressions,  the  superscripts  dm  and  am 
(derote  (quantities  associated  with  waves  of  dilatation  and  waves 
of  distortion,  respectively,.  The  definitions  of  q^  and  q^  are 
those  appropriate  to  a  viscoelastic  medium  as  given  in  Eq.  ,(3-4-16)  . 
In  the  case  of  .an  elastic  medium,  these  expressions  reduce  to 

eu  ■  bk  ^  and  n,  ■  bk  .  where  k  „  and  k  .  are  (defined  in 
3  P3  4  P4  P3  p4 

the  text  immediately  following  >Eqs.  (3.3.1)  and  (3.3.3),  respec¬ 
tively. 


5-  Superposition  of  Solutions  to  Enforce  the  Boundary  Conditions 

The  procedure  in  this  section  parallels  that  followed  in 
Section  III. 8-  We  combine  the  contributions  presented  in  the  pre¬ 
ceding  sections  so  as  to  satisfy  the  boundary  conditions,  namely, 
that  the  inner  boundary  of  the  cylinder  be  traction  free,  and  that, 
at  the  outer  boundary,  the  cylinder  displacements  and  appropriate 
stress  components  ate  continuous  with  those  of  the  surrounding 
medium.  Thus, 

at  p  »  P  (inner  boundary  of  the  cylinder). 
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Td_  +  t8_  -  0  ,  (4.5.1b) 
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and  at  p  *  1  (outer  boundary  of  the  cylinder). 
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where,  again-:  1)  the  superscripts  4  and  s  are  affixed  to 
-quantities  stemming  from  the  dllatational  and  -shear  -solutions 
to  the  wave  equation  in  the  cylinder,;  :2)  the  .simper scripts  dm 
and  sro  are  appended  to  stresses  .and  displacements  resulting 
from  the  dilatational  and  shear  -solutions  to  the  .wave  equation 
in  the  medium;  and  3)  the  superscript  8  identifies  the  -boundary 
tractions  and  -displacements  associated  with  the  incident  wave. 


After  the  terns  entering  Into  <4. 5.1)  are  computed  from  the 
appropriate  expressions,  and  the  common  trigonometric  and  exponen¬ 
tial  factors  are  divided  out,  the  resulting  nonhotoogeneous  sys¬ 
tem  of  linear  equations  in  the  .unlcnown  coefficients  .of  super¬ 
position,  Z. _ ,  may  be  written  in  matrix  form  as 

j.,pn 


'ij,-pn 


{*i,-pn} 


•(4.5.2) 


(p  “  2,  3,  (n  m  0,  1,  2,  ...)  , 


Where  C- ,  .  _  stems  from  the  coefficient  of  (.  in  the  i 

i j $ P«  j  > pn 

equation  of  Eq.  (4.5.1);  pn  stems  from  the  term  on  the 

right-hand  side  of  the  equation  of  Eq.  (4.5.1);  ( £,  }, 

J  *  Pn 

the  solution  vector,  is  a  6x1  column  matrix,  whose  complex 
elements  are  the  unknown  coefficients  of  superposition.  It 
should  be  noted  again  that  there  is  a  different  matrix  equation 
(4.5.2)  for  each  pair  of  values  of  p  and  n. 
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Ihe  e  lements  of  C,,  .  and  3 .  are  shown  in  detail  in 

ij.pn  i>pn 

Appendices  A  and  B. 


Aft ci  SWXViLl^  J-UXT  {  S> 


j»pn 


)  for  a  sufficient  number  of  values 


of  p  and  n  to  ensure  adequate  convergence,  we  may  calculate, 
to  the  desired  accuracy,  the  stresses  and  displacements  in  the 
cylinder  as 


d  .  s  d  .  s 

T  •  T  +T  .T  ■  T  +1  . 

rr  rr  rr  ^  rg  r§ 


-  x?_  +  Tf.  ,  (4.5.3a) 
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d  ,  s 
v  -  v  +  v 
r  r  r 


d  ,  s 

v.  -  v.  +  v_  , 


(4.5.3b) 


where  the  terms  appearing  on  the  right-hand  side  of  the  above 
expressions  are  given  in  Eqs.  (4.3.1)  and  (4.3.2). 


SECTION  V 


SPECIAL  FEATURES  OF  THE  PRESENT  ANALYSIS 


l  .  Initial  Con H -I  ft  nn a 


We  recall  that  tha  solution  constructed  in  the  preceding  sec- 
tions  is  the  response  of  an  elastic  cylinder  to  a  periodic  train 
of  pulses.  However,  by  measuring  time  from  the  moment  of  pulse 
arrival,  we  may  obtain  the  response  to  a  single  pulse  acting  on 
a  cylinder  initially  unstrained  and  at  rest,  provided  that  we 
adjust  the  time  between  pulses  so  that,  to  a  suTHcient  degree 
of  accuracy,  we  achieve  conditions  of  vanishing  stresses  and  dis¬ 
placements  in  the  cylinder  and  its  vicinity  —  prior  to  the 
arrival  of  each  pulse.  Indeed,  the  success  of  the  present  tech¬ 
nique  depends  chiefly  on  the  possibility  of  obtaining  stresses 
and  displacements,  before  the  pulse  arrives,  which  are  negligibly 
small  compared  to  significant  stress  and  displacement  levels  at 
later  times.  Essentially,  this  end  is  attained  by  providing  suf¬ 
ficient  "rest  time,"  that  is,  making  the  time  interval  between 
pulses  long  enough  to  permit  the  cylinder  wall  to  radiate  suffi¬ 
cient  energy  out  into  the  medium  and  return  to  an  undeformed 
state  of  rest.  Because  the  half  period,  T,  is  arbitrary,  the 
parameters  that  control  the  rest  time  are  to  some  extent  arbi¬ 
trary,  enabling  us  to  make  an  optimum  choice  of  the  rest  time. 

Some  discussion  of  this  procedure  is  in  order  here.  The  rest 
time  parameter,  tg,  can  be  nondimensionalized  as  (t^c 

that  is,  as  the  ratio  of  the  rest  time  to  the  time  required  by 
the  incident  wave  to  perform  a  transit  of  the  cylinder.  By 
referring  to  Fig.  2,  it  can  be  seen  that  the  rest  time  parameter 
can  be  expressed  in  terms  of  two  additional  nondlmensional  parame¬ 
ters  that  appear  in  the  analysis,  (cWttT/ 2b)  and  (c^mA/ 2b) ,  as 


cdmt0  cdmT  cdmA 

b  2b  '  2b  ‘ 


(5.1.1) 


The  pulse  duration  parameter,  (c^A/Zb),  while  not  shown  ex¬ 
plicitly  in  the  development  preceding  Eq.  (5.1.1) ,  enters  into 
the  Fourier  coefficient  ap  that  determines  the  wave  form.  As 

Eq.  (5.1.1)  indicates,  when  che  pulse  duration  is  prescribed,  the 
rest  time  parameter,  and  hence  the  rest  time  between  pulses,  may 
bo  varied  by  varying  the  parameter  (c^mT/2b) .  This  stems  from 

the  fact  that  the  half  time,  T,  enters  the  analysis  as  an  arbi¬ 
trary  parameter.  Thus,  tr  any  computations,  the  stresses  and 
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displacements  are  obtained  accurately  only  when  by  means  of  trial 

r'Amm  f  ^  nr>  o  \  ra  1  no  r\  £  ( r*  'P/OK\  an/4  Uam  a  a  a  u4  ■«* 

. r - '  - - —  N~dnr' — '  9  - - 

has  been  so  selected  as  to  reduce  the  initial  stresses  and  dis¬ 
placements  of  the  response  to  a  sufficiently  low  level.  For¬ 
tunately,  it  has  been  found  in  the  many  cases  presented  that 
satisfactory  results  may  be  obtained  if  Cc<jmT/2b)  is  kept  in 

the  reasonably  narrow  range  10  <  (c^T^b)  <  100.  The  cases 

considered  covered  a  wide  range  of  cylinder-medium  impedance 
mismatch  and  cylinder  thicknesses,  as  well  as  a  variety  of  pulse 
shapes . 

Improved  accuracy  may  be  obtained  by  treating  separately  the 
response  for  small  and  for  large  values  of  the  nondimensional 
time  variable,  {c^EJlb),  where  f  -  (t  -  tg)  and  measures 

time  with  respect  to  the  arrival  of  the  incident  wave.  From 
Eq.  (5.1.1),  it  can  be  anticipated  that  a  decrease  in  the  parame¬ 
ter  (cdmA/2b)  would  lead  to  further  improvement  in  the  initial 

conditions.  Ihus,  for  times  less  than  the  pulse  duration,  it  may 
be  advantageous  to  assign  a  value  to  \c^./2b)  that  is  less  than 

the  specified  pulse  duration.  The  corresponding  computations  are 
valid  for  times  less  than  this  altered  pulse  duration  (for  later 
times  the  actual  pulse  duration  parameter  must  be  used) .  This 
technique  is  particularly  suitable  for  the  case  of  a  rectangular 
pulse,  since  in  that  case  no  change  in  the  value  of  ap  is  re¬ 
quired.  The  foregoing  considerations  provide  the  basis  for  the 
device  by  which  an  incoming  step  pulse  may  be  treated,  namely, 
by  making  the  pulse  duration  long  enough  so  that  the  response 
reaches  a  steady  state  while  cdm't/2b  remains  less  than  (c(JmA/2b)  . 

One  further  consideration  is  necessary  to  satisfactorily  treat 
the  initial  conditions  requirement.  Since  we  have  not  fixed  the 
cylinder  with  respect  to  any  frame  of  reference,  it  will  trans¬ 
late  in  response  to  each  of  the  periodic  pulses.  In  order  that 
the  solution  be  applicable  to  the  case  of  a  single  transient 
pulse,  we  must  remove  from  the  displacements  the  rigid  body  con¬ 
tribution  resulting  from  the  preceding  pulse.  That  is,  we  re¬ 
quire  that  before  the  pulse  arrives. 


ur  cos  6  -  u0  sin  9  *  0 


(5.1.2) 
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Referring  to  Fig.  2  we  note  that  we  may  apply  this  condition  at 
any  time  t  where  -tQ  <  t  <  tQ,  provided  that  the  interval 

(  1  t*nn^*aofinf  0  a  » 4  wia  IV  a  mamma**  J  -  .  .L.  J  _  1_ 

\  -Q#  —  /  *“  — I"  —  w-— —  “  — -  -  V  A..A.  •  4LIIV  »•—**»»«#*.  nu^vu 

the  requirement,  (5.1.2),  is  enforced  is  now  described. 


Using  the  dilatational  and  distortional  potentials  given  by 
(3.2.7)  and  (3.2.14),  we  may  write  the  displacements  in  the 
cylinder  by  using  (3.2.4)  and  (3.2.11).  Ihe  resulting  expres¬ 
sions  may  be  put  in  the  form 
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(5.1.3b) 


where  the  upn  and  upn  are  functions  of  the  radial  coordinate 
r.  Enforcing  (5.1.2)  with  the  displacements  (5.1.3)  evaluated 
at  t  -  t  requires  that 
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as  well  as 
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~  Y  Z  U0R  sin  nS  ®  P  *  0.(5. 1.5) 
p-1  n-2 


For  convenience,  the  quantities  uP^  and  uP^,  appearing 

in  Eq.  (5.1.4),  are  evaluated  at  r  -  a.  Hius,  after  substituting 
(5.1.4)  into  (5.1.3),  we  may  calculate  the  displacements  in  the 
cylinder  as 
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(5.1.6b) 


For  the  case  o£  an  Incident  shear  wave  we  may  cast  the  dis¬ 
placement  expressions,  in  a  form  similar  to  those  given  by 
(5.1.3),  as 
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(5.1.7b) 


+  L  cPi  ' 

p-1 

where  v£n  and  v|n  are  functions  of  the  radial  coordinate  r. 

T  u 

The  term  ^  °pi  has  been  Included  in  the  expression  for  v^  to 
account  for  the  possibility  of  a  rigid  body  rotation. 
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In  this  case,  since  the  particle  motion  is  perpendicular  to 
Che  direction  of  wave  propagation,  the  consequence  of  not  fixing 
the  cylinder  with  respect  to  any  frame  of  reference,  is  a  trans¬ 
lation  of  the  cylinder  in  response  to  each  of  the  periodic 
pulses  —  in  a  direction  normal  to  that  of  the  progression  of 
the  wave.  In  addition,  a  rotation  of  the  cylinder  may  occur. 

To  remove  the  translation  and  rotation,  we  require  that  the 
following  two  conditions  be  satisfied  separately: 

vr  sin  6  +  vg-  cos  e  -  0  ,  (5.1.8a) 


v5-0 


(5.1.8b) 


Again,  referring  to  Fig.  2,  we  see  that  we  may  apply  these  con¬ 
ditions  ’at  any  time  t  where  -tn  <  t*  <  tn,  provided  that 

/v*  u  *"  —  u 

the  interval  ("t0,  t  )  represents  a  sufficient  rest  time. 

Thus,  enforcing  (5.1.8),  with  the  displacements  evaluated  at 

_  . 

t  ■  t  and  r  -  a,  leads  to 
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2.  Convergence  Criteria 


Criteria  to  control,  the  degree  of  convergence  of  the  various 
summations  invoivea  in  calculating  the  stress  and  displacement 
quantities  were  Introduced  into  the  computer  program  used  in  the 
present  study.  These  criteria  are  new  described. 


The  stresses  and  displacements  may  be  represented  in  the  form 

00  OP 

Ojj  ■  y  (p*  t)  ,  ifj  •  tpQ  , 

p-1  n»0 

OO  00 

ul  m  Y  T  e>  fc)  »  i  -  r ,Q  . 

p«l  n«0 

it  it 

Integers  N  ,  N  ,  M  ,  relating  to  the  degree  of  accuracy 

'ft 

desired,  were  selected  and  the  summations  were  carried  out  so 
that: 

1)  For  each  integer  p#  and  for  N  consecutive  terms 
In  the  n  summation, 


.  Pk 
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,  r-,  pn 


,  pk 
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■w1 - <10  • 

^  pn 
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The  value  N  of  the  integer  n,  at  which  the  condition  required 
by  step  1)  is  met,  of  course  varied,  with  the  integer  p,  the 
point  (p,  ©,  t),  and  the  subscript  ij,  or  i. 
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2)  Using  the  appropriate  value  of  the  summation  over  n 
up  to  n  *  N,  the  summation  over  p  was  carried 
out  so  that  for  M*  consecutive  terms. 


p-1  N 

I  l  l  <$(P.e,t)| 


k*l  n*l 


N 


!  £  I  ^(p^pt)! 


k-l  n-1 


<  10 


Control  of  the  computations  was  maintained  to  evaluate  the 
influence  of  the  summations  over  n,  on  the  current  summation 
over  p>  so  that  the  n-summation  (number  of  circumferential 
modes)  could  be  either  terminated  or  increased  if  required. 

Codes  present  in  the  computer  output  indicated  whether  the  con¬ 
vergence  criteria  were  met.  All  the  results  presented  were 
obtained  with  at  least  graph  plotting  accuracy.  In  the  neighbor¬ 
hood  of  maximum  values,  accuracy  to  four  figures  was  obtained 
in  most  cases* 
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SECTION  VI 


DISCI1SSIQM  Of  NUMB&ICAI  &ESULTS 


l.  fisncrai  JEgnalaeaticaa 


An  IBM  digital  computer  was  employed  to  obtain  the  solution 
of  the  system  of  matrix  equations  (3.8.2)  and  (4.3.2)  and  to 
carry  out  the  associated  stress  and  displacement  computations 
defined  by  Eqs.  (3.8.3)  and  (4.5.3).  The  computer  program  is 
described  in  Volume  II  of  the  present  report  (Ref.  16).  In  the 
discussion  that  follows,  even  though  results  are  presented  for 
specific  dimensional  values  of  geometric  and  physical  parameters, 
these  results  are  applicable  for  the  wider  class  of  cases  covered 
by  the  nondimens ional  values  of  these  parameters,  namely,  the 

h  A  ya 

ratios  ?,  “7,  7™,  and  77.  This  applies  for  the  case  of  a 

*  ^  F  F*  7 

viscoelastic  medium  as  well  as  for  an  elastic  medium. 


It  should  be  pointed  out  that  the  results  presented  through* 
out  this  report  are  given  in  nondimens  ional  form.  For  example, 
in  Fig.  5  we  have  given  the  hoop  stress  and  displacement  time 


histories  in  terms  of  the  quantities 


and  — 

V  '  ao  b 


respec¬ 


tively.  The  sign  of  the  actual  stress,  or  displacement,  is  thus 
determined  by  the  sign  of  Oq,  the  stress  associated  with  the 


incident  wave.  In  the  case  of  the  response  to  an  incident  shear 
wave,  it  is  the  sign  of  Tq  that  determines  the  sign  of  the 


actual  stress  or  displacement. 


2.  Computations  Performed  to  Illustrate  the  Validity  of  the 
Present  Method 

In  this  section,  we  illustrate  the  capability  of  the  present 
technique  to  satisfy  the  requirements  of  an  initial  un deformed 
state  of  rest  in  the  cylinder.  For  this  purpose,  the  stress  and 
displacement  response  to  a  rectangular  pulse  lasting  five  transits 
of  the  cylinder  has  been  computed.  The  Fourier  coefficients  a  , 
which  define  the  incident  wave  form,  are  ” 
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(6.1.1) 


where  a 
P 


is  the  Lanczos  factor  (Ref.  17)  that  can  be  included 


as  a  factor  in  the  Fourier  coefficients  to  improve  convergence 
of  the  series  representing  the  incident  wave  in  the  neighborhood 
of  discontinuities.  When  a  is  set  equal  to  unity,  Eq.  (6.1.1) 


leads  to  the  usual  Fourier  series  representation. 


Ihe  results  of  the  computations,  and  pertinent  parameters 
employed,  are  shown  in  Fig.  5.  After  the  pulse  has  passed,  and 
as  the  disturbed  cylinder  radiates  energy  outward  into  the  medium, 
the  maximum  hoop  stress,  at  Q  -  90°,  is  seen  to  vanish. 

This  shows  that,  although  steady-state  solutions  are  employed  by 
the  present  technique,  an  unstrained  state  of  rest  can  be  attained 
in  the  cylinder  and  adjacent  medium  before  the  arrival  of  the  next 
pulse.  Note,  for  exampl  *,  the  circled  points  in  Fig.  5a.  these 
represent  the  computed  values  of  aQQ  in  the  Intervals 

-1  <  (c^mt/2b)  <  0  and  11.5  <  (c^TLJlb) .  Referring  to  Fig.  5b, 

we  note  that  during  early  transit  times  the  character  of  the  dis¬ 
placement  response  is  nonlinear  and  suggests  considerable  rapid 
deformation  of  the  cylinder,  later,  the  displacement  time  his¬ 
tories  of  the  representative  points  selected  become  linear  and 
parallel  —  suggesting  motion  essentially  of  rigid  body  charac¬ 
ter  —  and  in  response  to  an  infinite  step  pulse  would  continue 
to  exhibit  this  behavior.  After  the  rectangular  pulse  has  passed, 
the  cylinder  regains  its  original  shape,  but  remains  displaced  in 
the  amount  of  the  rigid  body  motion.  It  was  precisely  this  rigid 
body  displacement  that  was  removed  from  the  periodic  solution  in 
Section  V.l,  entitled  "initial  Conditions.4'  Thus,  Figs.  5a  and 
5b  illustrate  the  feasibility  of  the  present  technique. 


3.  Comparisons  with  Published  Data  for  the  Response  to  an 
Incident  Dilatatlonal  Step  Pulse 

Comparisons  of  results  obtained  for  limiting  cases  of  the 
present  problem  with  corresponding  published  data  are  exhibited 
in  Figs.  6  through  9.  In  all  comparison  cases,  the  incident 
wave  is  a  step  pulse.  As  stated  earlier,  an  incident  rectangular 
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pulse,  of  duration  sufficient  for  the  response  to  reach  static 
conditions,  can  be  substituted  for  an  incident  step  pulse,  ftiere- 
fore,  the  Fourier  coefficient  ap  has  the  form  given  in  (6.1.1). 

Figures  6  and  7  are  concerned  with  the  response  of  an  infinite 
cylindrical  cavity  in  an  elastic  medium.  In  Fig.  6,  the  cavity 
displacements  obtained  by  the  present  technique  at  9  -  0°,  90°, 
and  180°,  are  compared  with  those  given  in  Ref.  3.  Agreement  is 
excellent  except  for  early  times. 

The  maximum  cavity  stress  response  to  a  step  pulse  is  the 
hoop  stress  at  9  •  90°.  In  Fig.  7a,  values  of  this  stress  cal¬ 
culated  by  the  present  analysis  are  compared  with  those  obtained 
in  Refs.  2,  4,  and  5.  For  times  greater  than  the  first  half 
transit,  essential  agreement  prevails  between  the  prerent  method 
and  Refs.  4  and  5.  Hie  stresses  predicted  in  Ref.  2  are  somewhat 
lower  than  those  given  by  the  other  analyses  up  to  the  time  when 
maximum  values  are  reached.  Then  all  results  approach  the  static 
value  given  by  the  Kirs  eh  formulas  <Re£.  18).  The  results  of 
Refs.  4  and  5,  shown  on  Fig.  7a,  are  for  various  truncations  of 
the  representation  of  the  deformation  in  circumferential  modes. 

As  pointed  out  previously,  the  present  analysis  differs  from 
those  of  the  other  references  in  that  it  does  not  constrain 
a  priori  the  number  of  circumferential  modes  taken  to  represent 
the  solution.  Instead,  it  permits  this  choice  to  be  governed, 
at  each  space  and  time  point,  by  the  condition  that  the  solution 
converge  to  a  specified  degree  of  accuracy.  For  times  less  than 
the  first  half  transit,  there  are  significant  differences  between 
the  results  obtained  in  this  fashion  and  those  of  the  other  referen¬ 
ces.  By  contrast,  results  from  the  present  analysis  for  a  three¬ 
mode  truncation  <and  which,  therefore,  represent  less  accurate 
predictions  of  the  response)  are  Identical  with  the  corresponding 
values  obtained  for  the  same  truncation  in  Ref.  4  over  the  entire 
range.  Thus,  it  may  be  inferred  from  Fig.  7a  that,  for  the  case 
of  an  incident  step  pulse,  a  three-mode  representation  gives 
reasonably  satisfactory  values  during  times  in  the  neighborhood 
of  the  maximum  response,  but  constitutes  a  highly  inaccurate 
approximation  to  the  deformation  during  early  transit  time.  This 
is  shown  even  more  forcibly  in  Fig.  7b.  Nevertheless,  the  close 
agreement  that,  for  a  three-mode  truncation,  exists  between  the 
results  of  the  present  analysis  and  those  of  Ref.  4,  indicates 
that  the  time  dependence  is  treated  properly  by  the  present  tech¬ 
nique  . 
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In  Fig.  7b  the  hoop  stress  confuted  by  the  present  method 
at  9  -  0°,  for  the  case  of  a  step  pulse  impinging  on  a  cylin¬ 
drical  cavity,  is  compared  with  corresponding  data  in 

Refs.  2,  4,  and  5.  Again,  these  references  present  results 
obtained  by  various  truncations  of  the  representation  in  cir¬ 
cumferential  modes,  and  again,  some  degree  of  agreement  prevails 
among  these  computations  during  later  transit  times.  However, 
wide  disagreement  exists  for  earlier  times,  particularly  during 
the  first  transit,  when  even  the  sign  of  the  stress  response  is 
in  contention. 

As  in  the  comparison  at  Q  -  90°,  the  results  obtained  at 
0-0°  by  the  present  method,  for  a  three -mole  truncation,  agree 
very  closely  with  those  of  Ref.  4  over  the  entire  range.  Over 
most  of  the  range  considered,  and  especially  for  the  early  transit 
times,  these  results  are  at.  variance  with  those  obtained  when 
the  number  of  modes  is  adjusted  to  comply  with  the  requirements 
of  convergence  at  each  space  point  and  time  (being  considered. 

This  shows  that  the  discrepancy  stems  from  the  circumferential 
mode  representation  end  not  from  the  treatment  of  the  time  de¬ 
pendence.  It  Is  unlikely  that  a  deformation,  and  hence  a  stress 
response,  which,  'during  the  early  stages  of  transit  of  the  inci¬ 
dent  wave,  involves  only  part  of  the  cavity  {namely,  its  front), 
could  be  properly  represented  by  three  modes..  In  view  of  this, 
and  since  the  initial  value  predicted  by  the  present  method  for 
the  hoop  stress  at  0  *  0°  is  zero  {pointed  out  in  Ref.  4  to  be 
the  correct  value),  it  is  felt  that  the  early  response  is  more 
realistically  described  by  the  present  technique,  even  as  com¬ 
pared  to  Ref.  5,  In  Which  six  modes  were  employed.  After  one 
quarter  transit,  the  close  agreement  with  Ref.  5  is  to  be  noted. 


In  Fig.  8  the  results  of  the  present  method,  which  stem  from 
the  theory  of  elasticity,  are  compared  with  the  corresponding 
data  given  by  Ref.  7  fcr  the  hoop  stress  response  in  a  thin  cylin¬ 
drical  shell  of  various  thicknesses  (h/R  »  .0048,  .‘019,  .0381) 
at  0  ■  90°,  to  an  incoming  step  pulse.  Two  media  with  signifi¬ 
cantly  different  shear  modulii  are  considered.  The  results  of 
Ref.  7  stem  from  an  analysis  of  the  cylinder  by  a  thin  shell 
theory  in  which  the  cylinder  acts  as  a  membrane.  Fcr  each  shell 
curve,  two  sets  of  points  obtained  by  tn*-  present  technique  have 
been  plotted,  one  for  the  hoop  stress  at  the  middle  surface,  and 
a  second  at  the  inner  surface  representing  the  additional  effect 
of  bending.  As  can  he  seen,  the  difference  between  the  middle 
surface  and  inner  boundary  stresses  appears  to  be  negligible  for 
the  thinnest  of  the  three  shells.  For  the  thickest,  the  differ¬ 
ence  reaches  a  value  of  about  10  per  cent.  Moreover,  some 
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differences  exist  between  the  response  as  determined  by  the  present 
analysis  and  that  of  Ref.  7.  Of  course,  it  is  only  to  be  expected 
that  elasticity  theory  and  thin  shell  theories  will  yield  differing 
results.  However,  it  i s  observed  in  Fig.  8  that,  as  the  shell  be¬ 
comes  thinner,  the  shell  theory  predictions  do  not  come  into  closer 
agreement  with  those  of  elasticity  theory.  Therefore,  on  the  basis 
of  the  available  information,  it  19  not  possible  to  determine  at 
what  value  of  the  thickness  to  radius  ratio  the  thin  shell  theory 
ceases  to  be  valid.  In  an  effort  to  determine  the  source  of  the 
divergence  of  the  results,  as  the  shells  become  thinner,  between 
the  shell  theory  predictions  and  those  of  the  theory  of  elasticity, 
we  note  that  in  Ref.  7  the  deformation  of  the  shell  is  constrained 
to  three  modes.  In  addition  to  the  total  response.  Fig.  9  shows 
the  time  history  of  the  contribution  o£  each  mode  as  given  in 
Ref.  7  for  the  case  of  the  shell  with  thickness  to  radius  ratio 
of  .019.  Points  computed  by  the  present  method,  but  with  the 
same  truncation  as  in  Ref.  7,  for  the  total  response  and  its  com¬ 
ponent  modes  are  also  shown  in  Fig.  9,  along  with  the  total  re¬ 
sponse  (present  analysis)  that  satisfies  the  convergence  require¬ 
ments.  The  agreement  that,  except  for  early  transit  times,  pre¬ 
vails  between  the  two  responses  demonstrates  that  the  discrepancy 
between  the  present  results  and  the  case  of  Ref.  7  is  not  due  to 
an  insufficient  number  of  circumferential  modes  in  the  analysis 
of  Ref.  7.  Figure  9  shows  that  the  major  discrepancy  resides  in 
the  third  mode  (the  n  ■  2,  or  inextensional,  mode).  It  is 
further  noted  that  the  results  for  the  cavity  case  and  all  the 
or. her  shell  cases  considered  in  this  report  (including  some  of 
those  of  Ref.  7  shown  in  Fig.  8b)  indicate  that  the  number  of 
transit  times  required  to  approach  the  static  condition  does  not 
seem  to  be  significantly  dependent  on  the  material  as  well  as 
geometric  properties  of  the  shell  and  the  surrounding  medium. 

This  value  seems  to  remain  at  about  six  to  eight  transit  times. 

By  contrast,  the  response  reproduced  from  Ref.  7  in  Fig.  9  de¬ 
parts  from  this  general  behavior  as  the  shells  become  thinner. 

4  Results  Obtained  by  the  Present  Analysis  for  the  Response 
to  an  Incident  Dllatational  Wave  (Elastic  Medium) 


Figure  10  illustrates  the  influence  of  liner  thickness  on  the 
response  to  an  incoming  dilatational  «tep  pulse.  The  maximum 
stress  response,  the  hoop  stress  at  6  ■  90°,  is  shown.  When  the 
incident  wave  passes  from  the  medium  into  a  relatively  stiffer 
cylinder  (slow  granite  into  concrete),  the  effect  of  increasing 
liner  thickness  is  to  reduce  the  peak  liner  stress.  When  passage 
of  the  incident  wave  is  into  a  relatively  softer  cylinder  (fast 
granite  into  concrete),  increasing  the  liner  thickness  increases 
the  peak  stress.  The  appearance  of  a  more  oscillatory  response, 
with  increasing  thickness,  should  be  noted  for  the  softer  liner. 
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Ihe  influence  of  liner  thickness  on  Llw  displacement.  response, 
when  the  liner  material  is  stiffer  than  that  of  the  medium,  is 
illustrated  in  Figs  11a  and  lib  (where,  for  ease  of  comparison, 
the  curves  of  Fig.  5b  are  presented  again  in  Fig.  lib) .  An  in¬ 
coming  aiiatationai  rectangular  wave,  five  transit  times  in  dura¬ 
tion,  is  the  input  for  both  of  these  cases.  Ihe  liners  in  Figs.  11a 
and  11b  have  a  thickness  lo  radius  ratio,  n/k,  equal  to  .01  and 
.2,  respectively.  As  expected,  the  thicker,  and  hence  stiffer, 
shell  experiences  much  less  deformation  of  the  cross  section 
(measured  by  the  vertical  distance  between  the  two  u^.  curves)  . 

Also,  as  can  be  determined  from  Figs.  11a  and  11b  by  the  onset  of 
pure  rigid  body  motion  (parallel  portions  of  the  curves),  the  full 
deformation  is  realized  earlier  in  the  thicker  shell.  It  may  be 
noticed  that,  as  the  back  of  the  wave  passes  over  the  cylinder,  the 
front  of  the  cylinder  (0  *  0e)  experiences  the  sharper  variation 
in  displacement.  This  effect  is  accentuated  by  an  increase  in  the 
liner  thickness.  After  passage  of  the  wave,  the  thicker  3hell 
regains  its  original  shape  first.  As  expected,  the  total  transla¬ 
tion  of  the  shell  is  independent  of  the  thickness.  Additional 
results.  Fig.  11c,  suggest  that  when  the  liner  material  is  suffi¬ 
ciently  softer  than  that  of  the  medium,  a  change  in  liner  thickness 
does  not  appreciably  Influence  the  over-all  deformation. 

Figure  12  contains  results  obtained  by  the  presbnt  analysis 
that  apply  to  a  dilatational  step  pulse  traveling  through  ah 
Infinite  elastic  medium,  and  impinging  upon  an  elastic  cylinder, 
for  four  separate  sets  of  cylinder  and  medium  elastic  properties. 
These  were  selected  so  as  to  study  the  Influence  of  impedance  mis¬ 
match.  Curves  G-C  and  G-S  apply  to  an  elastic  environment 
having  the  properties  of  granite  and  lined  by  concrete  and  steel 
cylinders,  respectively.  In  case  G-C,  the  incident  wave  trans¬ 
mits  energy  from  a  relatively  "stiffer”  to  a  softer  elastic  en¬ 
vironment  (value  of  Young's  Modulus).  In  case  G-S,  the  dis¬ 
turbance  is  transmitted  from  a  relatively  softer  to  a  stiffer 
elastic  environment,  with  attendant  magnification  of  peak  stress 
in  the  cylinder.  Curves  Sa-C  and  Sa-S  apply  tc  an  elastic 
environment  having  the  properties  of  sandstone  and  lined  also  by 
concrete  and  steel  cylinders,  respectively.  Again,  the  magnifi¬ 
cation  of  the  peak  stress  is  greater  as  the  disturbance  passes 
from  a  relatively  softer  to  a  stiffer  medium.  In  this  case,  the 
difference  in  raagnif ication  of  peak  response  is  even  more  striking 
but  not  surprising,  since  sandstone  is  much  softer  than  granite 
and  has  practically  the  same  properties  as  concrete.  Hie  depen¬ 
dence  of  the  magnification  of  the  incident  stress  amplitude  upon 
impeoance  mismatch  has  been  demonstrated  previously  in  Ref.  7 
for  thin  shells.  In  this  connection,  the  impedance  ratio  parame¬ 
ter,  (y  c^/yCd)  *  ^as  keen  indicated,  for  the  case  corresp -’ud- 

ing  to  each  curve,  on  Fig.  12.  It  will  be  observed  that  decreas¬ 
ing  values  of  this  parameter  are  associated  with  increasing  values 
of  peak  response.  72 


The  effect  of  tht  stress  wave  form  upon  the  maximum  stress 
for  a  liner  embedded  in  an  elastic  medium  is  illustrated  in 
Fig.  13.  The  results  for  three  cases  of  incoming  dilatatlonal 
plane  waves  are  presented.  The  wave  forms  (se«  2b)  ere: 

1)  a  step  pulse;  2)  a  triangular  wave;  and  3)  a  linear  rise** 
exponential  decay  wave  form,  having  the  same  rise  time  as  the 
triangular  wave,  The  expressions  for  ap,  which  describe  wave 

forms  2)  and  3),  are  given  in  Appendix  B.  All  three  waves 
have  the  peak  value  Oq.  As  can  be  seen,  the  step  pulse  produces 

the  most  critical  response,  eventually  reaching  a  static  value 
of  2.37.  The  response  in  cases  2)  and  3)  are  identical  during 
the  rise  time.  The  value  of  the  hoop  stress  in  these  cases  con* 
tinues  to  increase  during  the  passage  of  the  peak  of  the  incident 
wave  across  the  liner,  and  decreases  after  the  peak  has  passed 
the  back  of  the  liner.  The  response  to  the  triangular  wave  does 
not  reach  nearly  as  high  a  peak  vame  and  ^'’raases  more  rapidly. 
This  effect  stems  from  the  rate  of  decay  parameter  that  ueflnes 
the  exponential  decay  and  for  this  case  gives  a  higher  value  of 
incident  stress  to  that  wave  form  than  that  of  the  triangular 
wave  at  corresponding  transit  times.  All  responses  should  give 
a  value  of  aQQ  equal  to  zero  at  less  than  1/2  transit.  The 

slight  discrepancy  from  this,  shown  by  the  response  to  the  step 
pulse  in  Fig.  13,  reflects  the  fact  that  it  is  more  difficult  to 
represent  a  step  function  by  a  Fourier  series  than  a  linear  rise. 


5.  Results  Obtained  by  the  Present  Analysis  for  the  Cylinder 
Response  to  an  Incident  Dilatatipnal  Plane  Wave  in  a 
Vlacgfllaatlc,  Medium 

Results  of  a  computational  effort  concerned  with  the  response 
of  an  elastic  cylinder,  embedded  in  a  viscoelastic  medium,  to  an 
incident  dilatatlonal  plane  wave,  are  presented  in  Figs.  14 
through  19.  The  values  of  parameters  defining  th'  are 

identical  with  those  of  the  concrete  cylinder  er*’  »i 

elastic  medium  discussed  in  the  preceding  subsr  enac¬ 
tion  with  Fig.  12.  The  viscoelastic  propertie:  lum  are 

described  by  the  standard  linear  solid  dlscussev  -on  III. 5 

(see  Fig.  4).  The  values  of  the  instantaneous  ex  parameters 

and  the  mass  of  the  medium  were  held  constant.  T  properties 
are  identical  with  those  of  sandstone  considered  /tously  as 
an  elastic  medium  in  the  cases  treated  in  Figs.  and  13.  They 
are  given  by  the  values  of  Em,  v  ,  and  ym  shown  in  the  figures 

along  with  the  corresponding  values  for  the  elastic  liner.  Also 
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kept  constant  was  the  value  of  the  viscoelastic  parameter. 


1^/ ^2* 


which  was  set  equal  to  unity.  The  three  remaining  nondltnenslonal 
ratios  that  describe  the  viscoelastic  properties  of  the  medium  «ere 
permitted  to  range  through  a  set  of  values.  The  parameters 

and  *2/^2*  equal  to  each  other  in  each  case,  took  on  the 

magnitudes  .25  and  .1.  Ihe  relaxation-time  parameter, 

was  assigned  the  spectrum  of  values  0.  .1,  1,  10,  and  ».  Note 
that  the  lower  end  of  this  spectrum,  c<jm/^2^>  "  corresponds 


to  the  limiting  case  of  an  elastic-relaxed  material,  while  the 
upper  end-  c<jn/^2t>  *“  00 *  aPpi^es  to  another  limiting  case,  that 

of  an  elastic-unrelaxed  material. 


It  is  well  known  that  the  model  used  in  the  present  analysis 
describes  only  the  general  features  of  the  actual  viscoelastic 
behavior  of  soils  under  dynamic  loading.  Accordingly,  the  values 
assigned  to  the  viscoelastic  parameters  should  not  be  construed 
as  being  fully  representative  of  specific  materials.  However, 
these  values  are  significant  in  that  they  fall  within  practical 
ranges  and  lead  to  a  spectrum  of  results  from  which  meaningful 
qualitative  conclusions  can  be  drawn. 

The  influence  on  the  response  of  the  cylinder  of  variation 
in  the  viscoelastic  parameters  of  the  medium,  when  the  time  history 
of  the  Incident  stress  wave  at  the  front  of  the  cylinder  is  a  step 
function,  is  illustrated  in  Figs.  14  and  15.  The  maximum  stress 
response,  at  the  inner  surface  and  0  «  90°,  is  plotted  in 

Fig.  14.  In  these  figures,  the  upper  and  lower  curves  correspond 
to  elastic  media  whose  properties  are  described  by  the  unrelaxed 
Ccdm/^2^3  "  and  re^axed  (Cdn/*V>  "  ^  elastic  modulii, 

respectively. 

In  general,  the  effect  of  viscoelasticity  in  the  surrounding 
medium  is  to  increase,  in  some  cases  drastically,  the  stress 
response  in  the  liner  over  that  associated  with  the  elastic- 
unrelaxed  medium.  This  increase  is  characterized  by  two  major 
trends.  One  trend  is  associated  with  a  decrease  in  the  relaxation 
time  (indicated  by  a  decrease  in  the  value  of  <Wn2b>  >  the 

other  with  a  decrease  in  the  ratios  and  ^2/^2'  ^  early 

times,  before  the  maximum  response  associated  with  an  elastic- 
unrelaxed  medium  is  reached,  the  increase  in  stress  corresponding 
to  a  decrease  in  the  relaxation  time  is  relatively  slight  for 
larger  values  (of  order  100)  of  the  parameter  ^ut 
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becomes  far  more  significant  as  the  magnitude  of  the  relaxation 
time  decreases  to  the  order  of  one  transit  time  of  the  elastic*' 
unrelaxed  wave  front  (c^/fi^b  -  10,  1) .  This  can  be  seen, 

for  example,  in  Fig.  14a,  which  corresponds  to  the  case 
Tl^ftl  "  t2^2  "  •21’*  UP  to  the  time  when  the  peak  stress  is 

reached  (c{jmt/2b  <  2.5),  the  stress  response  in  the  case  c^/f^b-lOO 

is  identical  with  that  corresponding  to  a  purely  elastic  (unre~ 
laxed)  medium.  During  the  same  time  interval,  the  stress  for 
Cdn/^2b  “  ^  8rows  t0  a  value  about  20  per  cent  larger  than  the 

peak  stress  in  the  case  of  the  purely  elastic  medium  while,  for 
cdn/fi2b  "  reac^es  a  value  100  per  cent  larger.  As 

cdn/fi2b  the  8C7ess  approaches  the  value  associated  with  the 

purely  elastic -relaxed  medium,  which  is  more  than  double  that  of 
the  peak  stress  in  the  liner  embedded  In  the  purely  elastic- 
unrelaxed  medium.  It  should  be  noted  that  the  peak  stress  in 
the  case  of  the  elastic-relaxed  medium  is  almost  three  times 
that  in  the  case  of  the  elastic-unrelaxed  medium.  This  pattern 
of  increased  stress  response  for  decreased  relaxation  time,  ob¬ 
served  at  c(jmt/2b  -  2.5,  is  magnified  with  time  as  the  material 

in  the  viscoelastic  environment  continues  to  relax  and  the  liner 
must  resist  larger  deformations  in  the  medium.  Ultimately,  the 
stress  can  be  expected,  in  all  the  cases  where  the  medium  has  a 
finite  relaxation  time,  to  reach  the  steady-state  conditions 
associated  with  an  elastic-relaxed  medium.  This  process,  clearly 
Indicated  in  Fig. 14,  is  shown  to  be  extremely  slow  for  large 
values  of  <Wn2b’  100  <  cdn/^2b  ^  H°wever,  the  approach 

to  the  steady-state  condition  occurs  far  more  rapidly  as  the 
value  of  c^/^b  decreases  below  10.  When  this  condition  is 

reached,  the  stress  has  about  2-1/2  times  the  value  of  the  peak 
stress  associated  with  the  elastic-unrelaxed  medium.  Setting 

■  ^2/^2  “  1  in  the  creep  laws  developed  in  Section  III. 5 

is  another  way  to  reduce  these  laws  to  the  stress-strain  rela¬ 
tions  pertaining  to  the  medium  whose  properties  are  purely  elas¬ 
tic  (unrelaxed) .  Thus,  when  and  T2^2  aPProac^  unity, 

the  viscoelastic  curves  for  the  entire  spectrum  of  relaxation 
times  must  merge  into  a  common  curve  —  the  curve  associated  with 
the  elastic-unrelaxed  medium.  Taking  that  curve  as  a  basis  for 
comparison,  the  sizable  spread  of  the  viscoelastic  curves  in 
Fig.  14a  gives  a  clear  indication  of  the  magnitude  of  the  magni¬ 
fication  experienced  by  the  major  stress  in  the  liner  due  to  a 
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decrease  in  the  value  of  Tj/^  and  Fig111®  14b,  in  which 

Tl^fil  "  t2^2  “  *1'  gives  further  evidence  of  this  trend.  It  can 

be  noticed  in  this  figure  that  the  peak  stress  in  the  case  of  the 
elascicTuiaxeu  medium  is  now  five  times  greater  chan  che  peak 

stress  in  the  case  of  an  elastic -unrelaxed  medium,  as  compared 
to  three  times  greater  when  Tj/Oj.  -  T2/f>2  »  .25  (Fig.  14a). 

Figures  14a  and  14b  both  show  that  the  curves  corresponding  to 
the  two  elastic  media  diverge  quickly  from  each  other.  Since 
the  curves  corresponding  to  finite  relaxation  times  range  in 
between  these  two  curves,  it  may  be  concluded  that,  if  the  ratios 
of  and  T2/fl2  are  appreciably  less  than  unity,  the  repre- 

sentation  of  a  viscoelastic  medium  purely  by  its  elastic-unrelaxed 
or  its  elastic “relaxed  properties  may  lead  to  results  that  are 
either  unconservative  or  too  conservative. 

In  the  present  discussion,  valid  conclusions  can  also  be 
drawn  for  step  pulses  of  finite  duration.  Depending  on  the  pulse 
duration  and  the  degree  of  viscoelasticity  in  the  medium,  stresses 
may  be  Induced  in  the  burled  cylinder  that  are  drastically  larger 
than  the  peak  stress  predicted  by  considering  the  medium  to  be 
purely  elastic  (unrelaxed) .  The  longer  the  duration  of  the  pulse 
or  the  smaller  the  relaxation  time,  the  higher  the  stress  will 
be,  with  the  elastic-relaxed  condition  constituting  an  upper 
bound.  Conversely,  short  pulses  in  a  medium  with  a  sufficiently 
long  relaxation  time  may  not  induce  a  response  appreciably  differ¬ 
ent  from  the  response  associated  with  the  purely-unrelaxed  medium. 
As  a  final  note  concerning  Fig.  14  it  should  be  mentioned  that 
the  discrepancies  between  the  curves  at  the  time  of  arrival  of 
the  pulse  can  be  attributed  to  difficulties  sometimes  encountered 
in  approximating,  by  means  of  Fourier  series,  the  discontinuity 
of  the  Incident  step  pulse  for  the  values  of  the  parameter  b/c.T 
necessary  to  obtain  a  sufficient  rest  time  (see  Section  V) . 

Figure  15  is  concerned  with  the  deformation  of  the  inner 
surface  of  the  cylinder  (described  by  u^  at  0°,  90°,  and  180°) 
due  to  the  incident  step  pulse.  As  expected,  the  trends  shown  by 
the  curves  in  this  figure  parallel  those  for  the  stress  response. 

It  is  seen  that  the  principal  effect  of  viscoelasticity  of  the 
medium  is  to  Increase  the  deformations  over  those  predicted  for 
a  purely  elastic  (unrelaxed)  medium,  without  changing  their 
character  appreciably.  The  deformation  Increases  with  decreasing 
relaxation  times  or  decreasing  ratios  of  elastic-relaxed  to 
elastic-unrelaxed  properties  T2/n2).  The  full  liner 

deformation,  indicated  by  the  onset  of  rigid  body  motion,  is 
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reached  later  when  the  medium  is  viscoelastic?  than  when  the  medium 
is  purely  elastic  (relaxed  or  unrelaxed) .  The  rigid  body  motion 
itself  increases  when  the  relaxation  t<m«  nr  tv.®  ratios  of  relaxed 
to  unrelaxed  elastic  properties  are  reduced. 

The  influence  of  viscoelasticity  of  the  surrounding  medium 
on the  response  of  the  cylinder  to  an  incident  stress  wave  with 
a  triangular  time  history  at  the  front  of  the  cylinder  (see  Fig.  13) 
is  shown  in  Figs.  16  and  17.  During  the  period  of  increasing 
stress  response,  the  family  of  curves  in  Fig.  16  shows  a  magnifi¬ 
cation  of  the  stress  with  decreasing  relaxation  time  that  is  very 
similar  to  that  observed  in  the  case  of  the  incident  step  pulse. 

The  peak  stress  associated  with  the  elastic-relaxed  medium  is 
nearly  2-1/2  times  the  peak  stress  associated  with  the  elaatlc- 
unrelaxed  medium.  Again,  the  peak  stresses  corresponding  to 
finite  relaxation  times  in  the  medium  are  bracketed  between  these 
two  stresses.  Consequently,  results  obtained  by  treating  a  visco¬ 
elastic  medium  simply  as  an  elastic -unrelaxed  or  elastic-relaxed 
medium  may,  for  this  incident  wave  form,  again  be  either  unconserva¬ 
tive  or  too  conservative . 

It  is  Interesting  to  note  the  combined  effect  that  the  rise 
time  of  the  incident  wave  (1/2  transit)  and  the  value  of  the  re¬ 
laxation  time  in  the  medium  seem  to  have  in  the  response  of  the 
cylinder.  It  is  seen  In  Fig.  16  that  early  in  the  response, 
during  a  time  interval  comparable  to  the  rise  time  of  the  inci¬ 
dent  wave,  and  for  a  medium  With  relaxation  time  at  least  as  long 
as  the  rise  time,  the  stress  is  identical  with  that  associated  . 
with  the  elastic-unrelaxed  medium.  It  should  also  be  noted  that, 
after  the  passage  of  the  wave  over  the  cylinder  >4), 

the  liner  stress  decays  more  slowly  in  cases  involving  surrounding 
viscoelastic  media  than  in  cases  involving  elastic  media.  An 
explanation  for  the  behavior  of  the  stress  response  in  the  case 
of  a  viscoelastic  medium  is  that  it  must  involve  primarily  the 
slow  process  of  creep  strain  recovery  in  the  medium. 

The  deformation  of  the  liner  inner  surface  due  to  the  inci¬ 
dent  triangular  pulse  is  shown  in  Fig.  17,  and  is  again  described 
by  ux  at  0°,  90°,  and  180°.  For  the  sake  of  simplicity  we 

have  presented  only  the  curves  corresponding  to  the  elastic- 
relaxed  and  unrelaxed  media  and  to  the  viscoelastic  medium  for 
which  ^formation  f°r  fche  case  c^/^h  *  1.0, 

is  given  in  a  subsequent  figux^  (Fig.  19c).  It  is  seen  in  Fig.  17 
that,  as  in  the  case  of  an  incident  step  pulse,  the  effect  of 
viscoelasticity  in  the  medium  is  to  magnify  the  deformations 
without  changing  their  essential  character.  The  lower  the  relaxa¬ 
tion  time  the  larger  the  deformation,  with  the  elastic-relaxed 
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curves  constituting  an  upper  bound.  In  the  viscoelastic  case, 
after  the  wave  has  passed,  it  is  interesting  to  note  how  very 
long  it  takes  for  the  cylinder  to  regain  its  original  undeformed 
snape.  inis  is  inaicacea  oy  the  slow  converging  of  tne  curves 
after  the  maximum  deformation  has  been  reached. 

We  now  discuss  another  point  of  interest.  Referring  ahead 
to  Fig.  19c,  and  in  conjunction  with  Fig.  17,  we  note  the  attenu- 
ation  of  the  maximum  acceleration  with  decreasing  values  of  the 
relaxation  time  down  to  at  least  c^g/n2^  "  which  the  inner 

surface  at  Q  ■  0°  experiences  just  after  the  passage  of  the 
pulse  Cc<jm‘E/2t>  £  5).  However,  the  sharpness  of  the  peaking  of 

the  displacement  response  at  6  *  0°  for  the  case  of  the  elastic 
relaxed  medium  indicates  that,  at  some  value  of  the  relaxation 
time  in  the  range  0  <  c^m/ fl^b  <  1.0,  the  maximum  acceleration 

would  start  to  increase  again  and,  for  lower  values  of  the  re** 
laxat ion  time  down  to  c^/fijh  “  0,  would  exceed  that  corre¬ 
sponding  to  the  elastic-unrelaxed  medium.  v 

The  influence  of  incident  wave  history  on  the  response  of 
a  liner  embedded  in  a  viscoelastic  medium  is  illustrated  in 
Figs.  18  and  19..  Figure  18,  which  is  concerned  with  the  maximum 
stress  response,  combines  the  results  plotted  in  Fig.  14a  for 
the  case  of  an  Incident  step  pulse  with  those  plotted  in  Fig.  16 
for  the  case  of  a  pulse  of  triangular  form.  It  also  includes 
results  for  the  pulse  with  the  linear  rise-exponential  decay 
history  described  in  Fig.  13  and  Sedtion  VI. 4.  For  the  sake  of 
comparison,  the  curves  presented  in  Fig.  13  are  reproduced  in 
Fig.  18a  for  the  case  where  the  medium  is  purely  elastic  (un¬ 
relaxed)  .  The  only  other  figure  in  this  set  containing  results 
for  a  linear  rise-exponential  decay  pulse  is  Fig.  18c,  for  the 
case  c^j/fi2b  "  1.0. 

It  is  seen  from  these  figures  that  the  effect  of  the  visco¬ 
elasticity  in  the  medium  is  to  accentuate  the  differences  between 
the  stress  responses  to  the  three  pulses.  These  differences 
Increase  with  decreasing  relaxation  time.  It  can  also  be  seen 
from  the  curves  in  Figs.  18a  through  18d  that  the  difficulties 
encountered  in  computing  the  onset  of  the  stress  response  to  a 
step  pulse  were  not  encountered  in  the  case  of  the  other  two 
pulses . 

A  comparison  of  the  deformations  of  the  inner  boundary  of 
the  cylinder  due  to  the  three  incident  stress  waves  is  made  in 
Fig.  19.  This  figure  combines  the  curves  presented  in  Figs.  15 
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and  17.  It  also  presents  the  results  obtained  for  the  linear 

rise ‘exponential  decay  incident  pulse  In  the  case  of  the  elastic- 

unrelaxed  medium  and  the  medium  for  which  e.  /0-h  *  1.0.  v™- 

am  t. 

all  the  degrees  of  viscoelasticity  considered,  the  step  pulse 
induces  the  largest  cylinder  deformation  (as  measured  by  com¬ 
paring  the  vertical  distances  between  the  uf  curves)  and  tha 

most  rapid  translation  (indicated  by  Ug  at  90°).  The  triangu¬ 
lar  pulse  Induces  the  smallest  deformations.  Of  the  two  pulses 
of  finite  duration  (triangular  and  linear  rise-exponential  decay), 
the  triangular  pulse  causes  the  smaller  final  rigid  body  displace¬ 
ment  of  the  cylinder  and  that  displacement  is  attained  earlier. 
Conversely,  the  peak  acceleration  during  the  decay  of  these  two 
pulses  shown  to  occur  at  6  ■  0®  by  the  sharper  peaking  of  the 
corresponding  ur  curve,  is  the  greatest  in  the  case  of  the 
triangular  pulse.  Figure  19  demonstrates  that  the  viscoelasticity 
of  the  medium  has  the  effect  of  magnifying  the  differences  in  the 
response  to  the  various  pulses.  This  magnification  la  accentua¬ 
ted  by  reduction  of  the  relaxation  time. 


6.  Comparisons  with  Published  Data  of  the  Results  Obtained 
by  the  Present  Analysis  for  the  Response  to  an  Incident 
Shear  Wave 

In  Fig.  20,  the  hoop  stress  response  to  an  incident  shear 
wave  at  7  -  45°  and  135°  is  compared  with  the  results  pre¬ 
sented  in  Refs.  4  and  19.  The  remarks  made  in  the  discussion  of 
Fig.  7  apply  here  as  well.  At  times  corresponding  to  the  maximum 
stresses,  the  various  results  do  not  differ  too  widely,  and  all 
approach  the  static  value  given  by  the  Kirsch  formulas.  However, 
for  early  times,  when  the  incoming  shear  wave  is  still  in  transit 
across  the  cavity,  three  circumferential  modes  do  not  accurately 
describe  the  deformation.  This  again  can  be  seen  by  comparing 
with  the  other  references  the  results  obtained  by  the  present 
analysis,  where  we  used  that  number  of  modes  sufficient  to  ob¬ 
tain  convergence  of  the  solution  and  also  a  three-mode  truncation. 
When  the  present  computations  were  restricted  to  three  modes,  it 
can  be  seen  that  the  results  duplicated  those  of  Ref.  4.  However, 
the  more  accurate  solution  differs  from  that  of  Ref.  4  in  that 
the  maximum  is  reached  at  S’  -  135°  at  an  earlier  time  than  at 
e  -  45°. 
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7 .  Results  Obtained  by  the  present  Method  for  the  Response 
to  an  Incident  Shear  Wave  (Elastic  Medium) 

I.ie  cases  presented  m  Figs.  7.1  and  22  were  selected  so  as 
to  study  the  influence  of  liner  thickness  on  the  response  to  an 
incident  shear  wave.  The  liner  and  medium  parameters  are  the 
sane  as  those  selected  to  obtain  the  results  shown  in  Pig.  10  for 
the  case  of  an  incident  dilatational  wave.  The  .ia j or  stresses 
are  the  circumferential  stresses,  t_  at  6  «  45°  and  135°. 

ee 

As  was  noticed  in  the  cavity  comparison  case  (Fig.  20),  the  cir¬ 
cumferential  stress  at  &  »  135“  reaches  its  peak  more  rapidly 
than  the  corresponding  response  at  9  »  45°.  The  stresses  at 
both  angular  positions,  for  these  shells,  apparently  reach  the 
same  maximum.  For  these  cases,  the  response  per  unit  amplitude 
of  the  incident  wave,  as  well  as  the  overshoot,  is  greater  than 
the  corresponding  response  to  an  incident  dilatational  wave. 

These  effects  were  also  observed  in  the  cavity  response  (Refs.  4 
and  17) .  As  in  the  case  of  an  incident  dilatational  wave,  an  in¬ 
crease  of  liner  thickness  decreases  (or  increases)  the  response 
when  the  liner  is  stiffer  (or  softer)  than  the  medium. 

In  Fig.  22,  the  influence  of  liner  thickness  on  the  displace¬ 
ment  response  la  shown.  Effects  similar  to  those  noticed  in  the 
response  to  an  incident  dilatational  wave  are  indicated.  As  ex¬ 
pected,  the  path  of  the  rigid  body  motion  is  perpendicular  to 
the  direction  of  motion  of  the  incident  wave  front.  Also,  as  in 
the  case  of  the  incident  dilatational  wave,  when  the  back  of  the 
wave  passes  over  the  cylinder,  the  front  of  the  cylinder  (8  -  0°) 
experiences  the  sharper  variation  in  displacement. 

In  Fig.  23,  the  influence  of  imp e dance  mismatch  on  the  re¬ 
sponse  to  an  incident  shear  wave  is  illustrated.  The  parameters 
of  the  liner  end  surrounding  medium  are  Identical  with  those 
utilized  to  obtain  the  results  presented  on  Fig.  12.  The  major 
stress  response,  t__,  occurs  at  8  -  45s  and  135°,  the  maxi¬ 
mum  at  0  ■  135°  being  reached  earlier  in  time.  It  was  noted  in 
discussing  Fig.  10,  and  is  again  observed  here,  that  the  magnifi¬ 
cation  of  peak  stress  is  greater  when  the  wave  travels  from  a 
softer  into  a  relatively  stiffer  medium.  Also  to  be  observed  are 
the  indications  that  the  response  to  the  shear  wave  involves  con¬ 
siderably  greater  deformation  effects  than  the  corresponding 
response  to  a  dilatational  wave.  For  example,  note  the  sharpness 
of  the  peak  response  on  curve  Sa-S,  and  also  the  tendency  for 
the  stress  to  continue  to  oscillate  afterward. 
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The  influence  of  wave  form  on  the  response  to  an  incident 
shear  wave  is  shown  in  Fig.  24.  The  liner  and  medium  parameters 
are  the  same  as  those  for  WMf»h  rh*  rocnit; «  12  were  eb~ 

tained  for  the  incident  dilatational  wave.  The  major  stress 

response  is  that  of  the  hoop  stress,  t _ ,  and  occurs  at 

ee 

|  -  45°  and  135°,  with  the  maximum  being  reached  earlier  at 
S  -  135°.  The  general  character  of  the  response  is  similar  to 
that  for  the  case  of  an  incident  dilatational  wave.  The  most 
severe  response  again  occurs  when  the  incident  wave  form  is  that 
of  a  step  pulse;  and,  again,  for  the  parameters  considered,  the 
linear  rise-exponential  decay  wave  form  induces  a  hoop  stress 
response  of  greater  magnitude,  for  comparable  times  —  after  the 
passage  of  the  peak  of  the  incident  wave  —  than  does  a  triangu¬ 
lar  wave  form.  The  steeper  character  of  the  sheaf  wave  response 
has  already  been  mentioned  in  connection  with  previous  figures. 

The  abrupt  change  in  curvature  of  the  response  curve  at  9  -  135 % . 
for  the  case  of  the  step  pulse,  before  the  maximum  is  reached 
should  be  noted.  The  curves  at  9  -  135“  corresponding  to  the 
other  wave  histories  exhibit  a  similar  but  more  pronounced  trend. 
In  particular,  the  curve  associated  with  the  linear  rise- 
exponential  decay  pulse  shows  a  well-defined  double  peak. 


8.  Results  Obtained  by  the  Present  Analysis  for  the  Cylinder 
Response  to  an  Incident  Shear  Wave  in  a  Viscoelastic  Medium 

The  numerical  results  for  the  response  of  the  elastic  cylin¬ 
der  to  an  incident  shear  wave,  where  the  cylinder  is  embedded  in 
a  viscoelastic  medium,  are  presented  in  Figs.  25  through  30.  They 
were  obtained  for  the  same  cylinder  and  viscoelastic  media  as  in 
the  case  of  the  incident  dilatational  wave  {Section  VI. 5).  The 
histories  of  the  incident  shear  waves  are  also  the  same  as  for  the 
dilatational  waves.  They  are  described  in  Fig.  24.  To  assess  the 
influence  of  the  viscoelasticity  of  the  medium,  results  for  the 
elastic-unrelaxed  medium  are  included  as  a  limiting  case.  These 
results  have  already  been  presented  in  Figs.  23  and  24,  and  dis¬ 
cussed  in  the  preceding  section. 

The  major  stress  response  in  the  liner,  t_  at  the  inner 

ee 

surface  and  9  ■  45°  and  135°,  is  shown  in  Fig.  25  for  the  case 
of  the  incident  step  pulse  and  In  Fig.  26  for  the  case  of  the  tri¬ 
angular  pulse.  In  general,  these  results  indicate  that  the  prin¬ 
cipal  effect  of  viscoelasticity  in  the  surrounding  medium  is  again 
to  magnify  the  stress  response  as  compared  to  that  associated  with 
the  elastic-unrelaxed  medium.  The  trends  at  9  ■  45c  —  increasing 
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stress  response  with  decreasing  relaxation  time  and  also  with 
ratios  of  relaxed  to  unrelaxed  elastic  properties  —  are  very 
similar  to  those  at  9  -  90°  relating  to  the  incident  dila- 
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taining  to  Figs.  14  and  16  applies  equally  to  the  curves  for 
6  m  45°  in  Figs.  25  and  26a.  Tn  Flgi  25a  *  the  per  cent  in* 
crease  in  stress  at  c^t^b  -  2.5  resulting  from  the  decrease 

in  the  relaxation  time  of  the  medium  is  practically  identical 
with  that  observed  previously  in  Fig.  14a  at  c^t^b  -  2.5. 

Again,  the  ratio  of  the  peak  stress  associated  with  the  elastic- 
relaxed  medium  to  that  associated  with  the  elastic -unrelaxed 
Increases  from  a  value  of  about  3  in  the  case  “  t  2/fl2  -  .25 

(Fig.  25a),  to  a  value  of  about  5  in  the  case  t^/Gj  ■  .10 

(Fig.  25b).  The  hoop  stress  response  at  0  »  135°,  however,  de¬ 
parts  significantly  from  that  at  9  »  45®  in  the  case  of  both 
an  incident  step  pulse  and  an  incident  triangular  pulse.  In  both 
these  cases,  a  complete  reversal  of  the  trends  of  stress  magni¬ 
fication  with,  viscoelasticity  in  the  occurs  at  ~9  «  135° 

during  a  time  interval  after  the  arrival  of  the  pulse.  In  the 
cases  \faexe  t,/G-  -  t./G.  »  .23  <Figs.  23a  and  26a) ,  this  de¬ 
crease  in  the  stress  response  with  decreasing  relaxation  time 
stops  before  the  stress  associated  with  the  elastic -unrelaxed 
medium  has  reached  its  maximum  {c£Bt/2b  -  2)  and  gives  way  to 

the  more  prevalent  trend  of  increasing  stress  response.  However, 
it  is  seen  in  Fig.  25b,  wtwre  «  x2^2  *  ,i'  that  t*1® 

hoop  stress  at  9  »  135°  in  the  case* of  the  elastic-relaxed 
medium  remains  considerably  below  the  peak  stress  for  the  elastic- 
unrelaxed  medium  until  about  a  half  transit  past  this  petit.  This 
indicates  that,  at  least  for  the  case  of  an  incident  step  pulse, 
the  effect  of  decreasing  the  ratios  and  is  not 

only  to  magnify  the  trend  itself  but  also  Its  duration.  As 
cdTr/G?b  a  peculiarity  associated  with  this  trend  at 

0  »  135°  is  the  appearance  in  Figs.  25a  and  26a  of  a  stress  of 
opposite  sign  for  a  short  time  after  the  arrival  of  the  pulse. 

Ihis  stress,  barely  noticeable  in  Figs.  25a  and  26a,  becomes  much 
more  prominent  as  t^/G^  and  ^2^2  decrease  (Fig.  25b/. 

Other  effects  of  the  viscoelasticity  of  the  medium  on  the 
stress  response  at  9  -  135°  can  be  seen  in  Figs.  25  and  26. 
flie  curve  representing  the  response  to  an  incident  step  pulse  in 
the  purely  elastic  (unrelaxed)  medium  goes  through  abrupt  changes 
in  curvature  before  reaching  its  peak.  For  the  case  of  a  triangu¬ 
lar  pulse  (Fig.  26),  these  changes  in  curvature  are  more  pronounced. 
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and  the  curve  has  what  amounts  to  a  double  peak.  With  a  decrease 
in  relaxation  time,  this  behavior  is  progressively  less  in  evi¬ 
dence.  It  ceases  to  be  noticeable  for  cHrn/fi7'o  <1.0. 

It  was  also  pointed  out  in  the  preceding  section  that  the 
stress  response  to  an  incident  step  pulse  In  the  elastic-uurelaxed 
medium  reaches  its  maximum  earlier  at  6  m  135°  than  at  9m  45°. 
It  can  be  seen  in  Fig.  25  that  this  is  also  the  case  when  the 
elastic-relaxed  medium  is  involved.  By  contrast,  the  hoop  stress 
response  associated  with  the  viscoelastic  media  for  which 
cdn/^2b  "  *-s  never  lar8er  at  9  m  135“  than  at  9  •  45°. 

The  influence  of  the  shear_wave  history  on  the  stress  response 

of  the  liner  inner  surface  at  6  m  45°  and  135°  is  illustrated 

in  Fig.  27  for  the  case  of  the  medium  for  which  c .  /fi-b  =  1.0. 

°  urn  2 

This  influence,  when  the  medium  is  viscoelastic,  can  be  compared 
to  that  when  the  medium  la  purely  elastic  (unrelaxed)  by  refer¬ 
ring  back  to  Fig.  24.  For  the  finite  value  of  the  relaxation 
time  considered  here,  it  can  be  seen  that  the  nature  of  the  dif¬ 
ferences  between  the  responses  to  the  three  distinct  pulses  de¬ 
scribed  in  Section  VI. 7  is  preserved.  Again,  the  effect  of 
viscoelasticity  is  to  magnify  these  differences.  It  can  also  be 
seen_that  the  disappearance  of  the  double  peak  for  the  response 
at  dm  135°  due  to  the  viscoelasticity  in  the  medium  occurs  in 
the  case  of  the  linear  rise-exponential  decay  incident  wave  as 
well  as  in  the  case  of  triangular  incident  wave.  By  contrast  with 
the  responses  shown  in  Fig.  24  for  the  case  of  the  purely  elastic 
(unrelaxed)  medium,  the  stress  responses  at  9  »  45°  in  the  case 
of  the  viscoelastic  medium  are  almost  identical  in  character  with 
those  at  9  ■  135°. 

The  nondimensional  displacement  response  perpendicular  to 
the  path  of  the  incident  wave,  at  9  -  0°,  90° 

and  180°,  and  at  the  liner  inner  surface,  is  presented  in 

Figs.  28  through  30.  It  is  given  by  v_  at  0°  and  180°  and 

9  — 

vf  at  90°.  Since  vr  is  an  odd  function  of  9 ,  the  distance 

between  opposite  points  along  the  diameter  perpendicular  £o  the 
path  of  the  incident  wave  remains  unchanged  and  vr  at  6  «  90°  — 

the  middle  curve  in  each  set  of  three  given  in  Figs.  28  through 
30  —  represents  the  translation  of  the  whole  cylinder.  When  the 
vertical  distances  between  this  curve  and  tut.  other  two  curves 
(v_  at  0°  and  180°)  are  equal  —  as  they  are  for  a  substantial 

9 

portion  of  the  response  —  the  diameter  of  the  cylinder  aligned 


with  the  path  of  the  wave  experiences  a  pure  rotation  in  addition 
to  a  translation  equal  to  that  of  the  perpendicular  diameter. 

The  curves  in  Figs.  28  through  30  are  very  similar  in  charac¬ 
ter  to  those  shown  in  Figs,  id,  i/,  and  iy  tor  the  response, 
2imx/oQb,  to  incident  dilatational  pulses.  They  indicate  similar 

tnnds.  For  all  three  incident  waves  —  the  step  pulse,  the  tri¬ 
angular  pulse,  and  linear  rise-exponential  decay  pulse  —  the 
effect  of  a  decrease  in  the  relaxation  time  is  to  increase  the 
translation  described  by  vr  at  6  ■  90°  and  the  rotation  de¬ 
scribed  by  the  other  two  displacements.  Figure  28  indicates  that, 
at  least  for  the  case  of  an  incident  step  pulse,  the  same  effect 
is  induced  by  a  decrease  in  the  ratios  t^/£^  and  T2^2* 

curves  in  Fig.  29  for  the  case  of  an  incident  triangular  pulse 
indicate  that  the  sharp  acceleration  of  the  front  of  the  cylin¬ 
der,  just  before  the  intensity  of  the  incident  pulse  at  this 
point  in  the  cylinder  reduces  to  zero  (ctmt/2b  -3),  is  at 

first  attenuated  by  decreasing  values  of  the  relaxation  time.  Fig 
ure  30c  shows  that  this  trend  continues  at  least  until  cdn/fi2b 

reduces  to  the  value  1.0.  However,  as  implied  by  the  character 
of  the  curve  for  the  displacement  at  9  ■  0°  in  the  case  of  the 
elastic -relaxed  niedium,  at  some  value  of  c^/J^b  below  unity 

the  acceleration  starts  to  increase  again.  It  also  becomes 
sizable  at  ?  ■  90°.  The  very  slow  recovery  of  the  liner's  origi¬ 
nal  shape  after  the  passage  of  the  pulse  is  again  Indicated  in 
Fig.  29  for  the  case  by  the  slow  convergence  of 

the  curves. 

Figure  30  illustrates  the  influence  of  the  Incident  wave 
history  on  the  displacement  response.  Again,  decrease  in  re¬ 
laxation  time  magnifies  the  differences  between  the  responses  to 
the  various  pulses.  It  should  be  noted  that  for  the  case  of  a 
viscoelastic  medium  (Fig.  30c)  as  well  as  for  the  ct-e  of  the 
elastic-unre taxed  medium  (Fig.  30a),  the  sharp  acceleration  ex¬ 
perienced  by  the  front  of  the  liner  in  the  case  of  an  incident 
triangular  wave  is  far  less  pronounced  in  the  case  of  the  linear 
rise-exponential  decay  incident  wave. 
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.  Results  Obtained  by  the  Present  Analysis  for  the  Superposition 
of  the  Cylinder  Responses  to  Incident  Dilatational  and  Shear 
Waves  In  an  Elastic  Medium 

As  formulated  in  Section  I,  the  idealized  problem  with  which 
the  present  analysis  is  concerned,  is  that  of  determining  the 
response  in  a  cylindrical  liner  buried  in  an  infinite  medium,  the 
cylinder  being  subject,  in  succession,  to  plane  dilatational  and 
shear  waves.  Results  obtained  by  the  present  analysis  are  valid 
for  different  physical  situations,  provided  that  effects  other 
than  those  postulated  can  be  neglected.  For  example,  an  air- 
induced  ground  loading  results  when  a  pressure  wave,  generated  by 
an  above-ground  burst,  travels  along  the  surface  of  an  elastic 
half  space  with  superseismic  velocity,  V,"  where  V  >  c^.  Then, 

the  results  of  Ref.  20  imply  that  plane  dilatational  and  shear 
waves  are  transmitted  into  the  half  space.  The  transmitted  wave 
fronts  are  inclined  to  the  surface  by  the  angles  ■  arcsin  (c|n/V) 

and  at  -  arcsin  (ctm/V)  respectively.  These  fronts  are  inclined 

to  each  other  by  the  angle  9  (see  Fig.  1)  where  4>  ■  -  a^. 

Thus,  the  incident  shear  wave  Impinges  upon  the  liner  at  an  angle 
9  with  respect  to  the  dilatational  wave,  and  its  time  of  arrival 
is  delayed  by  an  amount  depending  upon  the  location  of  the  cylin¬ 
der  with  respect  to  the  surface  loading.  Now,  if  the  liner  is 
located  at  a  depth  such  that  the  effects  upon  its  response,  of 
surface  waves,  as  well  as  reflections  from  the  surface,  can  be 
neglected,  then  the  application  of  the  present  idealized  analysis 
is  Justified. 

Another  physical  situation  for  which  the  present  analysis 
could  be  applicable  was  suggested  in  Ref.  21  and  involves  the 
direct-induced  loading  case  (burst  in  direct  contact  with  surface) . 

In  this  case,  justification  for  the  use  of  the  present  analysis 
requires  the  validity  of  certain  approximations.  These  Involve 
the  location  of  the  cylinder  with  respect  to  the  burst  such  that, 
in  the  ground  effects,  a  spherical  dilatational  wave  front  may  be 
replaced  by  a  plane  dilatational  wave  front  (the  previous  assump¬ 
tions  on  surface  waves  and  reflections  are  also  required) .  If 
the  half  3pace  is  assumed  to  be  homogeneous  and  isotropic,  only 
the  response  to  an  incident  dilatational  wave  is  needed  to  approxi¬ 
mate  this  case.  If  the  liner  is  assumed  to  be  located  in  the 
lower  level  of  a  two-layered  half  space,  then  dilatational  and 
shear  waves  will  be  refracted  through  the  interface  (estumed  to 
be  plane  in  this  case).  The  angle  between  refracted  dilatational 
and  shear  waves  may  be  calculated  in  terms  of  the  incident  angle 
and  the  seismic  velocities  of  the  two  layers.  The  delay  between 
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arrive  1  times,  at  the  cylinder,  of  the  refracted  dilatational  and 
shear  waves  may  be  calculated  from  the  known  location  of  the  cylin- 
der  with  respect  to  the  layered  medium  interface.  Sufficient  dis- 
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may  be  neglected,  is  required. 

Two  numerical  problems,  posed  in  Ref.  21  to  illustrate  the 
conditions  discussed  above,  were  analysed.  The  parameters  involved 
have  little  significance  beyond  providing  somewhat  realistic  ex¬ 
amples  (on  empirical  grounds,  the  amplitude  of  the  shear  wave  was 
chosen  as  one-third  the  amplitude  of  the  dilatational  wave)  . 

As  an  example  of  air-induced  loading,  a  concrete  cylinder, 

50  ft  in  diameter,  is  regarded  as  embedded  in  a  homogeneous  slow 
granite  environment,  at  a  depth  that  requires  the  transmitted 
waves  to  travel  a  distance  of  150  ft  before  contact  with  the  cylin¬ 
der.  The  surface  wave  is  assumed  to  be  a  step  pulse,  with  velocity 
V  -  8500  ft/sec.  From  the  analysis  of  Ref.  20,  it  may  be  deduced 
that  Oq  and  xQ  are  negative  with  respect  to  the  cylinder  co¬ 
ordinate  systems  of  Fig.  1.  Results  obtained  for  this  case  are 
illustrated  irk  Fig.  31.  From  the  data  presented  on  these  figures, 
the  angle  qp  between  wave  fronts  is  19.7°,  and  the  time  of 
arrival  of  the  shear  wave,  with  respect  to  the  arrival  of  the 
dilatational  wave,  c^Bt/b,  is  calculated  as  c^Bt/b  -  4.38. 

On  these  figures  are  shown  the  components  of  the  hoop  stress 
response  at  the  inner  surface  to  the  dilatational  and  shear  waves, 
as  well  as  the  total  response.  Because  the  amplitude  of  the 
dilatational  wave  dominates,  the  maximum  response  for  the  points 
we  computed  is  obtained  at  9  ■  90°,  even  though  at  this  angle 
the  effect  of  the  shear  wave  response  is  to  reduce  the  stress. 

At  9  ■  64.7°,  corresponding  to  9  ■  45°,  the  shear  wave  re¬ 
sponse  also  reduces  the  stress,  while  at  9  m  154.7°  (©  ■  135°), 

an  addition  to  the  dilatational  response  occurs  after  the  arrival 
of  the  shear  wave. 

As  an  example  of  a  direct  induced  loading,  a  steel  cylinder, 

50  ft  in  diameter,  is  regarded  as  embedded  in  a  fast  granite 
layer,  separated  by  a  plane  interface  from  a  layer  having  a 
seismic  dilatational  velocity  c^.  “  6000  ft/sec.  A  plane  dila¬ 

tational  wave  (the  approximation  co  the  spherical  wave)  impinges 
on  the  interface  at  an  angle  of  incidence  qp^  ■  15°.  The  cylin¬ 
der  is  assumed  to  be  located  in  the  fast  granite  layer  at  a  depth 
such  that  the  refracted  waves  travel  100  ft  before  impinging  on 
it.  Denoting  the  refraction  angles  for  the  transmitted  dilata¬ 
tional  and  shear  waves  as  9^  and  the  angle  9  between  the 
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transmitted  wave  fronts  may  be  found  by  utilizing  the  date  in  Fig.  32, 
along  with  the  relations  sin  tP<j/c(im  “  8*-n  “  8*-n  <*)tm/ctm-  *n 

thi9  manner  the  angle  <p  is  found  to  be  <p  ■  22°.  The  delay  time 
of  the  refracted  shear  wave  with  respect  to  the  refracted  dilata- 
t tonal  wave  is  found  to  be  c^m6t/b  «  2.92.  Results  for  this  case 

are  presented  in  Fig.  32  in  which  component  dilatational  and  shear 
responses,  as  well  as  their  superposition,  are  shown.  It  should  be 
noted- that  the  dilatational  component  is  identical  with  that  shown 
in  Fig.  12.  The  shear  component  represents  that  shown  in  Fig.  23, 
with  a  different  abscissa  scale,  cdmt:^'  rather  than  the  corre¬ 
sponding  shear  WAve  parameter,  ctmt/b .  Again,  because  of  the 

dominance  of  the  dilatational  wave  amplitude,  the  maximum  hoop 
stress  for  the  points  we  computed  occurs  at  0  m  90°. 
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RECTANGULAR  WAVE  FORM 
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TRIANGULAR  WAVE  FORM 


LINEAR  RISE' EXPONENTIAL  DECAY  WAVE  FORM 
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Fig. 2b  RECTANGULAR,  TRIANGULAR,  AND  LINEAR  RISE- 

EXPONENTIAL  DECAY  WAVE  FORMS 
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ELEMENT  of  linen 
Showing  stpess 

NOTATION 


h°iri 


OILATATIONAL  WAVE 
PROGRESSING  WITH 
SEISMIC  VELOCITY 
cdm 


(a)  OILATATIONAL  WAVE 


ELEMENT  OF  LINEN 
SHOWING  STRESS 
N  OVATION 


To— H 

!  i| 


SHEAR  WAVE 
PROGRESSING  WITH 
SEISMIC  VELOCITY 
*t  m 


<b  )  SHEAR  WAVE 


P  i «  3  CYLINDRICAL  LINER  ,  INCOMING  WAVES ,  COORDINATE  SYSTEM, 
STRESS  AND  DISPLACEMENT  NOTATION. 
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Fig  4 
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MECHANICAL  MODEL  OF  A  STANDARD  LINEAR  SOLID 
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LINER  CONCRETE  E  •  2  3  *  10*  p»l ,  V  •  .  2  ,  y  •  A  .5  slug*  /  M.3 

MEDIUM  '  SLOW  GRANITE  ,  Fm  *  10  *  10*  p*i ,  |/m  ..2S  .  ym  .  3.2  »iu9»  /  f<  4 


O  Computed  point*  illutl'Olmg  tolilloclion  of  mitioi  condition*  cnttrio. 


< b>  DISPLACEMENTS 


Fig.  5  RESPONSE  OF  LINER  INNER  SURFACE  TO  INCIOENT 

OILATATIONAL  RECTANGULAR  WAVE  (  Eloilic  rr«dium) 
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0  4  0  Points  obtoined  by  present  onolysis. 


Fig.  6  COMPARISON  OF  CAVITY  DISPLACEMENTS  WITH 
THOSE  OBTAINED  BY  REF.  3. 

<  ELASTIC  MEDIUM) 


O  Results  obtained  by  present  method ,  using  number  of  modes 


required  by  convergence  criteria 
A  Results  obtained  by  present  m«thn»i>  «•»«**  r.*C,!,2 


«dmt/2° 

Note  Curves  reproduced  from  Ref.  4 

Fig.  7o.  COMPARISON  OF  PRESENT  COMPUTATIONS  WITH  THOSE  OF 

REFS.  2.4  ANO  3  FOR  HOOP  STRESS  AT  0  -90*  DUE  TO  INCIDENT 
OILATATIONAL  STEP  PULSE,  AT  BOUNDARY  OF  CAVITY  IN  ELASTIC 
MEDIUM  ,  •.  •  1/3 
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O  Results  obtained  by  present  method ,  using  number  of  modes 
required  by  convergence  criterio 

A  Results  obtained  by  present  method,  using  n»0,l,2 


Fig.  7b  COMPARISON  OF  PRESENT  COMPUTATIONS  WITH  THOSE  OF 

REFS  2  ,4,  ANO  5  FOR  HOOP  STRESS  AT  •  »  O*  DUE  TO  INCIDENT 
DlLATATlONAL  STEP  PULSE  ,  AT  BOUNDARY  OF  CAVITY  INELASTIC 
MEDIUM,  «  ■  1/3. 


LINER'  •£-«  019,  E  •  30  (10  6  )  psl ,  V  •  .25  y  •  1 52  slugs  / « t  3 
ME0IUM  ■  £m  ■  9(10*  )  psi ,  l/m  *  25  ,  ym  •  5. 2  slugs/ ♦  t.3 

□  A  Results  obtoined  by  present  method  ,  ot  0*  ond  90*  using  number  of  mode* 

r  «nn  Ku  «*mu*»**  »  *  • 

-/  - w,  .  •••(  •  V 

V  ©  Results  obtotned  by  present  method,  at  0*  ond  90*.  using  n*  0,1,2 


O  A  0  Results  obtained  by  present  method  for  n»  0, 1, ond  2. 


Note'  Curves  ond  asymptotes  reproduced  from  Ref.  7 
Fig  9  COMPARISON  OF  HOOP  STRESS  AND  ITS  MODAL  CONTRI¬ 
BUTIONS,  AT  LINER  MIDDLE  SURFACE,  WITH  THOSE  GIVEN 
IN  REF.  7. 
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Scole  applies  to  tost  granite 


LINER  ‘  CONCRETE  ,  E  *  2  5  *IO®psi,  v  »  2  ,  y *  4  5  slugs  /  «l  3 
MEDIUM  ■  SLOW  GRANITE  ,  Em  *  1.0  »  IO*pti ,  vm  *25,  ym  «  5  2  slugs  /fl. 3 


(b) 


Fig. II  DISPLACEMENT  RESPONSE  AT  LINER  INNER  SURFACE  TO 
INCIDENT  OILITATIONAL  RECTANGULAR  WAVE,  cdmA  , 

(Elastic  mtdiuml  ^ 
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LINES'  CONCRETE »  E  *2.5  *  <0*  p*i ,  V  .  2.  ,  y  •  4.5  itugt/Tt.5 
MEDIUM-FAST  GRANITE ,  Em  ■  9.0  a  10*  pal,  Pm  •  25,ym«5.2  aluga/tt.3 
-  '  -  —  "  -  ■■  ■  h  /  R  •  .0 1 


Fig.  lie  DISPLACEMENT  RESPONSE  AT  LINER  INNER  SURFACE  TO  INCIDENT 
DILATATIONAL  RECTANGULAR  WAVE  (Elastic  medium) 
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Fig.  12  INFLUENCE  OF  CYLINDER -MEDIUM  IMPEDANCE  MISMATCH  ON 
MAXIMUM  STRESS  IN  THE  CYLINDER  (0®$  ot  B  •  90°  ond  inner 
•  urfoeft)  FOR  AN  INCIDENT  DILATATIONAL  STEP  PULSE  IN  AN 
ELASTIC  MEDIUM.  ,  ns 


LINER :  CONCRETE  ;  h/R«.IO  t  E>  (2.5)110  Ipti  ;  V  «.20  ;  V  >4.S  slug*  /  ft 


F i<±  13  EFFECT  OF  INCIDENT  OILATATIONAL  WAVE  FORM  ON  MAXIMUM  STRESS  l  <T##  of  inner 
surface  and  e«90*l  OR  LINER  IN  ELASTIC  MEDIUM 
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.14  EFFECT  OF SURROUNOING  MEDIUM  VISCOELASTICITY  ON  MAXIMUM  STRESS  <  <T##  ol  90  *  end  Inn*' 

»u'foe«)  IN  LINER  FOR  INCIDENT  OILATATIONAL  STEP  PUISE 

10/ 


4.44  slugs /H 


INNER  SURFACE  TO  INCIOENT  DILATATiONAL  TRIANGULAR  WAVE 
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Response  to  triangular  save 

Response  to  lineor  rise  exponential  decay  wave 
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Fig.  19  EFFECT  OF  INClOENT  OILATATIONAL  WAVE  FORM  ON  DISPLACEMENT 
RESPONSE  AT  INNER  SURFAC^OF  LINER  IN  VISCOELASTIC  MCOIUM. 


It  r  g  n  o  rr.  b  r>  r 


O  Results  obtained  by  present  method 

modes  reduced  by  converger, it  criteria 


o  * 


A  Results  obtained  by  present  method,  using  n»  0,1, and  2 


Fig. 20  COMPARISON  OF  PRESENT  COMPUTATtONS  WITH  THOSE 
OF  REFS.  4  AND  18  FOR  HOOP  STRESS  (  ot  0  ■  45°  and 
135  •  )  DUE  TO  INCIDENT  SHEAR  STEP  PULSE,  AT 

BOUNDARY  OF  CAVITY  IN  ELASTIC  MEDIUM,  €  *  1/3 
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Scale  applies  to  fast  gra-tite 


LINER  CONCntU  ,  E  •  2  5ilOf'pi 
MEDIUM'  Slow  GRANITE,  Em  - 1 0« 


y  -  A  5  t  I'.Q!  /  f : 


i  IO®pti  .  «  25,  ym  «  5  2  s  l  u  9  *  /  M  5 


(  b) 

Fig.  22  EFFECT  OF  LINER  THICKNESS  ON  DISPLACEMENT  RESPONSE  OF  THE  LINER 
INNER  SURFACE  FOR  THE  CASE  OF  A  INCIDENT  SHEAR  RECTANGULAR  WAVE, 
(tloitie  matflum  I  elmA/26*7.5 
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Note.  Incident  wovi  form  given  in  Fig.  24 


fl.  25  EFFECT  OF  SURROUNDING  WEDIUM  VISCOELASTICITY  ON  HOOP  STRESS  AT 
INNER  SURFACE  OF  LINER  FOR  INCIDENT  SHEAR  STEP  PULSE 
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Note  ••  Incident  wav#  form  given  in  Fig.  24. 


Fig.  26  EFFECT  OF  SURROUNDING  MEDIUM  VISCOELASTICITY  ON 
HOOP  STRESS  AT  INNER  SURFACE  OF  LINER  FOR 
INCIDENT  TRIANGULAR  SHEAR  WAVE. 
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Fig.  27  EFFECT  OF  INCIDENT  SHEAR  WAVE  FORM  ON  HOOP  STRESS  AT  INNER 
SURFACE  OF  LINER  IN  VISCOELASTIC  MEDIUM  WITH  m,° 


Fig.  20  EFFECT  OF  SURROUNDING  MEDIUM  VISCOELASTICITY  ON  DISPLACE¬ 
MENT  RESPONSE  OF  LINER  INNER  SURFACE  TO  INCIDENT  SHEAR 
STEP  PULSE.  121 


LINER  =  CONCRETE  10,  E»(25KlO*)psi . 


LINER  INNER  SURFACE  TO  INCIDENT  TRIANGULAR  SHEAR  WAVE. 
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Not*  incidtnt  wo»*  form*  glvpn  In  Rig.  24 

pig  30  EFFECT  OF  INCIDENT  SHEAR  WAVE  FORM  ON  DISPLACEMENT  RESPONSE  AT  INNER  SURFACE 

OF  LINER  IN  VISCOELASTIC  MEDIUM. 
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CYUNOEft 


Fig  31  SUPERPOSITION  OF  CYLINDER  RESPONSE  (  Hoop  ttrooto*  ot  ii»«i«r  boundary)  TO  INCIDENT 
DILATATION AL  AND  SHEAR  WAVES  IN  AN  ELASTIC  MEDIUM. 


CYLINDER 


OIL  ATATlONAL  AND  SHEAR  WAVES  IN  AN  ELASTIC  MEDIUM. 
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B.  Representation  of  the  Fourier  Coefficients  Associated  with 
the  Incident  Wave 
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APPENDIX  A 


TRANSFORMATIONS  EMPLOYED  TO  FACILITATE  THE  rnMPtTTATtnMc 


Ir.  this  appendix  a  number  of  transformations  are  introduced 
to  facilitate  the  computational  effort  Involved  in  performing 
the  inversions  of  the  boundary  condition  matrix  equations  (3.8.2) 
and  (4.5.2).  For  ease  of  reference,  we  rewrite  these  equations 
omitting  the  Indices  of  summation  p  and  n,  but  understanding 
dependence  upon  them.  For  the  case  of  an  incident  dilatAtional 
wave,  we  have  Eq.  (3.8.2), 

h,]  {*,}  -  K} . 

and  for  the  case  of  an  Incident  shear  wave  we  have  Eq*  (4.4*2), 


Kd  {5j}  “  f  i] 


CA.2) 


The  stresses  and  displacements  associated  with  the  response 
to  the  incident  dilatations!  and  incident  shear  waves,  Eqs.  (3.8*3) 
and  (4.5.3),  respectively,  will  now  be  put  in  nondimensiooal  form. 
As  an  Intermediate  step  in  this  direction,  we  define 


■1J 


(A. 3a) 


’ij 


m  b 


i  •  r,  9  , 


(A.  3b) 


where  (a^,  u^)  are  presumes  to  be  calculated  by  Eq.  (3.8.3) 
and  (Tjj,  Vj)  by  Eq.  (4.5.3). 
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When  introduced  into  (A.l)  and  (A. 2),  the  transformation 
(A. 3)  nondiraensionalizes  the  matrix  elements  on  the  left-hand 
sides,  and  (A.l)  and  (A. 2)  may  be  written  in  the  form 


C 


..1 
“•j  j 


(A. 4) 


(A.  5) 


dilatational  and  shear  waves,  respectively,  and  where  C ^  and 


where  Cq  and  are  the  Amplitudes  of  stress  in  the  incident 

st 

are  notvH men « tonal.  Then,  if  we  let 

2 

bo* 

■  zj  -  (if)*)  - 


(A.6) 


The  boundary  conditions  become 


Kj3  ft}  -  {5t}  <  <*-7> 

and 

Pd  ft} '  ft}  •  <*•*> 


where  D,  and  ^  are  nondlmenslonal .  Referring  to  the  ex¬ 
pressions  in  Sections  Ill  and  IV  by  which  the  stresses  and  dis¬ 
placements  are  represented,  we  note  that  they  are  linearly 
homogeneous  in  the  Zj  and  Zy  Thus,  from  (A. 6),  when  the 

solution  vectors  to  (A. 7)  and  (A. 8)  are  introduced  into  the 
expressions  defined  by  (A. 3),  we  obtain  the  nondiraenisional 
stress  and  displacement 
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(A.  9a) 


2  n 


u . 


u. 


■f  rl  no  cn  ►U  / 


A  ~  -i  «4 - - 

*-»■*  V>  *-  c* 


--  -  1 

vi  4.  «vv*  wu  u  iv/ua  x. 


wave  ; 


awu 


-  „Iii 

«  To 


(A. 9b) 


in  the  case  of  the  incident  shear  wave. 


In  principle,  (A. 7),  (A. 8),  and  (A. 9)  may  be  used  directly 
to  determine  the  nondimen sional  stresses  and  displacements. 
However,  it  has  been  found  that  the  form  of  these  equations  may 
lead  to  computational  difficulties.  The  numerical  behavior  of 
the  Bessel  functions  J^(z) ,  Y^(z),  which  are  constituents  of 

c, , ,  c 


the  matrix  elements 

tlon.  For  fixed  z, 

small  end  the  Y_{z) 
n 


ij*  necessitates  a  further  modiflca- 

as  n  becomes  large,  the  Jn{z)  becomes 
becomes  large .  To  keep  -the  elements  of 

(A. 7)  end  (A.S)  within  the  bounds  of  the  IBM  7094  computer,  end 
to  ensure  calculations  of  a  well  behaved  character,  we  modify 
the  construction  of  Eqs.  (A. 7)  and  (A. 8).  To  emphasize  that  the 
elements  of  these  matrices  and  the  solution  vectors  depend  upon 
the  Indices  of  summation  p  and  n,  we  again  let  these  integers 
appear  explicitly.  Then  we  modify  the  elements  in  columns  1 
and  3  by  letting 


CU,pn  - 


2nn!  , 

Wkl“  lj'pn 


,n 


~  2  n*  — 

ij,pn 


id  Cil,„n  '(A'10a) 

'k1 


2j/Pn 


where 


(A.IOb) 


k  ■  1  when  J  ■  1  , 
k  ■  2  when  j  -  3  . 
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We  modify  tho  elements  in  columns  2,  A,  5,  and  6  by  lettinp 


kJn„ 


in 

q,.l 


ij>pn  2nnl  ij,pn  '  iJ>Pn  2nnl  • 


(A. 10c) 


•“ri...  .  C, 


1 


j,pn  |q  jn  j>pn  ’  j,pn  | q  |n  j,pn  * 
k  k 


(A.IOd) 


where 


k  «  1  when  j  “  2  , 

k  «  j  -  2  when  j  -  4 ,  5,6. 

After  the  definitions  (A.  10)  ere  imposed,  (A  -  7)  and  (A.  8)  become 


! 


9C 


ij 


N 


(A.  11) 


Kj]  {cj} "  f  4  • 


{A.  12) 


Hote  that  the  substitutions  defined  by  (A. 10)  have  modified  only 
the  left-hand  sides.  The  right-hand  sides  have  not  been  altered. 

Because  of  (A. 10),  instead  of  computing  the  Jn(s)  and 
YgCz),  In  columns  1  through  4,  we  compute  the  functions  defined 
by 


Jn(Z)  JnU)  '  Yn(2)-' 

*qk' 


2nni 


V*>  . 


(A. 13a) 


where  lc  ■  1  in  columub  1  auu  2;  and  k  ■  2  in  columns  3  and  4. 
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In  columns  5  and  6,  we  compute 


j  (2)  .  i-JU.  j  (z)  ,  Y  («)  -  Y(z) 

n  |z|n  n  n  2nnt  n 


(A.  13b) 


The  elements  composing  the  left-  and  right-hand  sides  of 
(A. 11)  and  (A. 12)  are  given  in  Appendices  B  and  C. 

For  computational  purposes,  Eqs.  (A. 11)  and  (A. 12),  which 
represent  sets  of  six  by  six  complex  equations,  were  cast  respec¬ 
tively  inu^  the  following  form,  which  represents  sets  of  twelve 
simultaneous  real  equations: 


PJ  ft}  ■  ft}  > 


(A.  14) 


ftj]  ft}  ‘  ft}  > 


(A.  15) 


where,  using  superscripts  R  and  I  to  denote  real  and  lmagi* 
nary  par^s,  (A. 14)  is  constructed  from  (A. 11)  as  follows: 


lu  - 


a  <  i  <  6),  a  <  s  <  6> 


LU  “  "  Cij 


(1  <  i  <  6),  (7  <  j  <  12) 


Li4«,j+6  "  Lij 


(1  <  i  <  6),  (1  <  j  <  6) 


Li+6,j  "  ‘  Lij 


(1  <  i  <  6),  (1  <  j  <  6) 


^  'Vrk 

V3? 


(1  <  J  <  6) 


'V  'VT 

Xj+6  “  Zj 


(1  <  j  <  6) 
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M1  -  D1 


-  Di 


(1  <  i  <T  6) 


(1  <  i  <  6)  • 


Similarly  (A. 15)  is  constructed  from  (A. 12)  in  the  following 
manner: 


U  ~1J 


"ij 


ij 


^i+6,j+6  "  ^ij 


i+6,j  "  ‘  ^ij 


_ 

‘  -i 


'i*.  -j 


"i+«  -  K 


(1  <  i  <  6),  (1  <  j  <  6) 

(1  <  i  <  6) ,  (7  <  j  <  12) 

(1  <  i  <  6),  (1  <  i  <  6) 

(1  <  i  <  6)>  (1  <  j  <  6) 

(1  <  j  <  6) 

(l  <  j  <  6) 

(1  1  i  <  6) 

(1  <  i  <  6)  . 


U 


We  recall  the  discussion  in  Section  IV. 1  and  note  that  L 

and  stem  from  the  response  to  Incident  dilp national  and 

shear  waves,  respectively.  Thus,  these  elements  may  fee  given  in 
the  same  form  if  we  regard  n  «  n  as  pertaining  to  L. ,  and 

_  ^  ^ 

n  «  -  n  as  pertaining  to  £ . , .  Then,  the  real  matrix  elements 

are  given  by 
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t-  ?  tr  e  c-  e  t+  e  c-  e 


L,1 


i)  -  -  ‘IPn^lW  -  35^5 


p) 


!  o 


-  -  q, 


[**)  -  «•  -  41w>  *  ^  WV> 

l‘VH"  <-(q  p)-  n  1  qp 


.  A. 

L 1 3  £ 


q2P 

(n-1)Jn(q2P)  “  2(2+1)  ^tv+l(q2P) 


L,4  a2 


■>  •% 

P 


(n-l)Yn(q2P)  -  2&(n+l)Yn+1(q2P) 


1,5  “  Ll,6  “  *•*  *  Ll,12  "  0 


T  m 

l2,l  p2 


q2P  o. 

2(n+l)  **trt-l^ql^  ‘  (n*1)Jn(qlP) 


L  •SL- 

L2,2  p2 


<\ <  q 

L 


2p(n+l)Yn+1(qiP)  -  (n-l)Yn(qiP) 
21 


-  -  21 
2,3  2 


{T^2~  ’  l)Vq2P)  +  FCn+T)  ^n+l(q2P) 


[_  (q2P) 


L  -  -  i 

L2,4  2 


(T^  '  1)  W>  ♦  ^  WV> 


(q,P) 


q2P 


r\,  f\* 


L2, 5  "  L2,6  "  *'•  “  L2,12  "  0 
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L3,l  '  '  ql 


L3,2  *  *  q! 


*©  ■  -  1} 
i©  -  Pr1  -  4 


VV  '  2(rin)  ^n+l<qi> 

Vql>  *  Y  i  (q  ) 

ql 


L3,3  "  n 


(n-l)Jn(q2)  “  2  (n+1)  ^n+l(q2J 


L3i4  "  n 


(n-l)Yn(q2)  -  2(n+l)?n+l(q2)j 


.n 


Oi  a ' 


<q3  -  q|> 


®fe-[«2>  +  (:DR3 


2(n+l)  q3 


(3^3  -  (§ 


*i 

j 


(n+l)3 


'  2 (Si)  q3[(^n+l)3  +  (j^Grtad 


+  2%q3 


2q3q3 


i© +  (§JV  B©  *  ©PS |  -  «<--«[(&%-({&] 


B© +  <- 


n3  12 


$  +  © 


~T 

Y. 


n3 


+  (i*  -  ?' 


)[*©♦©]%♦[$♦(?)& 


-(^H-^3^l^l)3  +  q3^U)3]} 
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t*  ? 


k 


L3,6  "  n 


(>?r  -  (X3 


rr 


■  '  (3?Si>)[vW  +  vUda!} 


,  nJl  ,  2in+1If-  «I  . :  > 

*  (n"1>Yn4  +  q4Y(n+l)4  q4r(n+l)4 


*4e 


rr 


(^7)  {<■»-«&  ■  (&lvU)*  •  %«*]} 


+  (”-1)^4  '  Vo 


(n+l)4  +  q4Y(n+l)4 


1 


L 


3,7  "  L3,8  "  L3,9  “  L3,10  “  0 


n3 


9,f(i)j(n+l)3  +  (1)J  (n+Ds] 


2(n+l)  H3 

+  2(n+l)  q3  [(S^n+l)  3  '  (f)J(rri-l)  3 1 

■  +  (Dfe  •  [2® +  ©R 

■  [(SR3 +  0R>]} 

-  (3  -  3)1  [2®  +  Qfc  -  [2CS)  +  (S)fi 
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*  -  (^)  {(^))[V(«H)J  '  VcrH-Ds]  -  tn"l)2n3 

^[^3  +  VU)3l-^ 


h.6  -  (H  •  (SK 


where 


Rr6  *  2  (^r)  '  <’4  -  q!>  Ftk  ’  2V4^4 


+  fr&r[V<,*.l)4  +  Vc*l)*]} 


+  [2n(n-l)  -  ttf “  q|)  ]yJ4  + 


v  lq4  <rH*1>4  V(t*fl)4 


■*Se  ‘  2  CSf.)  {I2"'"'1'  ‘  (’4  *  <4>  Pn4  +  2V4^4 


_  .4a..[~ 
(■rt-l)  lq' 


IvU.  ■  VU.I) 


-  |j(»Hn-l)  -  U*  -  ^ 


4(n+l)r- 


[q  y* 

l  V(i 


(n+1) 4  +  q4Y(n+l)4 


L4,7  “  L4,8  "  L4>9  “  L4,10  "  0 


137 


H,u  -  '  “[(SKs  +  (S4,] 


'4,12 


[i^v  +(h?  i 

LW  re  T  \m-/  re j 


n2 

q.  -v 


L5,l  '  nJn{,l)  "  2 (nil)  Jn+l(<1l> 


>•5,2  *  "W  '  2(“+1)Vl'’l) 


L,  _  ■  nJ  (a.) 
5,3  aNn2y 


L5,4  "  nYn*q2* 


5,5 


_n 


(■£■)  {°^3  *(2cJnr)[Vc»flj3+V(BH)j]}  +  <3 


2(0+1)  fr  Zl 


*3e 


[q3Y(n+l)3  "  q3Y(n+l)3 


5,6 


•‘"fe  & +  &} 


S,7  '  L5,8  ‘  L5,9  *  LS,10  '  0 
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r*  t 


J5,  11 


n 


(2a*-)  (, 

Vnl'  1 


nJn3  *  2 (n+1) /  q3J( 


(n+1)  3  '  VU,]}  -  <3 


5,12 


+  ziami  -  + :  ;i  s 

q3e  t^"*1)3  q3Y(tr+l)3j 

"fe  &  • 


~  jv 

L6,l  “  *  nJn(ql) 

L6,2  "  '  SW 

* . V 

6,3  "  ‘  “W  +  2 (n+1) 


_2 
q*>  ~ 


J6,4 


6,5 


“6,6 


"  «Yn(q2)  +  2(n+l)Yn+1(q2) 


-  (fe  R4  -  (5&))[vU4 + vU)} 


nYn4 

IV  (n+1)  4 

*v 

L6,8 

*v 

-W,9- 

**6,10  "  0 
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L6,12 


q.  r  ~j.  ~ 

V2(n+i))  i.<!4'T(n+i)4  (n+l)4 


!} 


-<* 


+  vU*] 


and 


•1J  "  ‘  Ll-6,  j« 


7  ^  i  i  12 
1  <  1  <  « 


LIJ  *  Ll-«,  J-* 


7*1*12 

7*1 


In  the  preceding  equations, 

*V 

Jnk-W 

•V 

\k  -  w 


k  -  3,  4  . 


The  expressions* in  real  form,  for  the  stresses  and  displace¬ 
ments  associated  with  the  responses  to  incident  dila^tional  and 
shear  waves  will  be  given  now.  The  transformations  leading  to 
Eq  (A  9).  in  conjunction  with  the  results  (A.IOo)  and  (A.  13). 
yield —  when  applied  to  the  appropriate  term.  In  the  exprea.ioe. 
defined  by  Eqs.  (3.8.3)  and  (4.5.3)  —  the  following  formulation 
of  the  stresses  and  displacements  associated  with  the  response  to 
an  incident  dilatational  wave.  Note  that  the  components  of  the 
solution  vector  entering  into  these  terms  are  those  obtained  from 
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Eq.  (A. 14).  In  addition,  the  rigid  body  contributions  discussed 
in  Section  V.l,  have  been  removed  by  taking  account  of  Eqs,  (5.1.6) 
and  (5.1.9),  and  the  variable  t/T  is  to  be  regarded  as  computed 
with  respect,  co  the  time  of  arrival  of  tne  incident  wave  as 


(!)  -  W  (*) +  (¥")  (^) 


Thus, 


-  -  cos  n©  |  q‘ 


(qrp) 


—  -  i  > 

2  V 


Pn(V>  - 


2(n+l)pJn+l<qlp)  r 


WV>  ^ 


-  i  <n-l)Jn(q2P)  “  2(n+l)  Jtrt-l(q2p)  X3 

.  P  ■  * 


-  K!) 

P  1 

. •  •  •  •»  • 

■  ql  ■  1)}Jn<,lf>)  *  2(n+l)p  Jn+l(qlp)X7 

+  4$©  ‘  (7772  *1)}^<v>  -  ^flf1  WV>  1*8 

L  (qjp)  J  q^p 


2 

q9p  ~ 


2  (n-l)Jn(q2P>  “  2(n+l)  Jirt-l*q2P^  X9 


■*2  (n-l)Yn(q2p)  -  2(n+l)  pYn+1(q2P)  ]xl0  sin  pir(*) 
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X  i 


~pn 

0Qd  -  -  COS  nd 


Cl /An  2  n(n-l)T  _ 
|2(^ql  +  2  jJn(qlp) 


+  _jL 


2(n+l) 


+  -^jw1  +  Wv> 


X, 


(n-l)Jn(q2p) 


-  -21 


q?p  - 


2<n+l)  rrt-1  vh2 


+  '  2(n+1)pW(V)R[CDS  K*). 


M  2  .  n(n-l)1'w 
ql  + 


]ja<V> 


ql  Wv> 


2(ttfl)  p 


*7 


^  +  +  WV>|** 


q2p 

+  ^  (n-l)Jn(q2p)  -  2{n+l)  Jn+l(q2p) 


+  ^[(n-l)Yn(q2P)  -  2<n+l)pYn+1<q2p) 


10 


sin  p*(£) 


:  I/" 
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pn 

Jr  ? 


sin  ri0 


JL 

p2 


q^p  ~ 

2  (u+1)  Jirt-l(qip' 


(n-l)  Jn  CqLP> 


X, 


2(n+l)pYn+L(qLp)  -  (n-l)Yn(q1p)]xj 


21 

hljininzll 

2  W)2 


1}VV>  +  to«37  WV> 


2» 

2 


-  2}w> 


l(v)* 


4(n+l) 

❖ 


Yn+l/q2p^ 


qlp  * 

2(n+l)  Jn+l(qlp) 


(n“l)Jn(qlP) 


+  -a- 

p2 


2(n+l)pYn+.(q  p)  -  (n-l)Yn(q1p) 


h 


W)2  ■i^,v 


(q  *p) 


(n+l) p  n+?  v^2 


2i 


(2n(n-l)  .  4 (n+l)  ~ 

{7~tt  •  +_xZY'>+i<q2p) 

L\q2p / 


V 


10 


cos  pir  £ 


sin  pv 

I 
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Hlrt 


'-pn 

u 

r 


+ 


+ 


+ 


+ 


+ 


where 


and 


COS  H@ 


{  p  Jn(qlp)  "  2 (nil)  Jn+l(qlp) 


.p  W>  ■  2^l)V1(V)]*2 


Lp 


Jn(q2P>]X3  V  [p  Vq2p)  k}jco^p7< t)  j 


{ 


_2 

q,  ~ 


I 


p  Jn<V>  *  2(^T>  WV>  V  [?  W>  '  2<“fl)  WV> 


[“■  w>  K +  [?  *«(y>'Ro}*taHi)  ] 


nl 


■{ 


J 2 

q-i  ~ 


p  Jl(ql«  -  4-  J2^1« 


xi  +  [p*w> 


[p 


x3  + 


{ 


J  W>  -  £  vv> 


J?1(q2P)]i4}co.[pXV)] 

X7  +  [p  W5  ’  «2(qlS) 


[p  VV>]x9 +  [p  V  V>  ]xio}sln(K 


* 

£_\  , 

T  /  J 


nl 


* 

T 


■{ 


0  if  n  ^  1 


if  n  -  1 
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~pn  rr  ^  1 A'  r""*  n'V 

u0  -  -  sin  ne  |[“  J^P)  jx2  +  [“  Yn  (q^lJXj 


+  p  Ji,(q2P>  “  2  (n+1)  J.i+l<q2p)  X3 

+  [?  Vv>  ‘  !(ltI)li¥  RHMJ)] 

/r**  <v  r™  n'V 

'  {[?  Jn<V>R  +  [?  W>]X8  +  ?  VV> 


_2 

q0  ~ 


2  (n+1)  Jn+l(q2p)  X9 


*  ip  Tn(V>  '  «»+l>WV>  RoHiKl)] 

ii  {[p  R<q&  R +  (i  w>R +  [i  W> 

+  [f  ^(qj?)  -  4Y2(q2B)jx4Jco«[PT  ^-j 


q2  1' 

4  W>J; 


.2 

q,  ~ 


ff  1  ~  "p  1  ~  1~  I  ~  4n  ~  ~ 

j[p  W>  K  +  [p  Yl(ql3)  JX8  +  8  Jl(q2«  ■  4  J2(’2« 


\  V1(q2P)  -  4Y2(q28)  jx10}sln[pi7 


where  6 


*  t 
and  ^  -  (-5 


|  0  if  n  ^  1 
l  1  if  n  -  1 

"  ^  '  2(w) 
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The  stresses  and  displacements  associated  with  the  response 
to  an  incident  shear  wave  are  obtained  from  the^preceding  expres¬ 
sions  bv:  l)  regarding  a.  u.  and  0  as  t.  v.  and  ■> . 

respectively;  2)  n  as  equal  to  -n;  3)  interchanging  sin  no 

*  /'“*  t  r 

_  t  t  0  -  / _ b  '  /  dm'.  , 

and  cos  no;  4)  replacing  —  by  “  -  "y  -  2  j;  and 

dm  tm 

5)  treating  the  n  »  0  case  separately  for  v^°  as 


{  ^  Vq2p)  X3  +  [21f1(‘l2P)lX4}‘:“s['”r  t! 


’2 

q0  ~ 


r  q,  ~  **  r  ~  i~  i  t 

-  >J  -f  JL(q2P>  X,  +  [2Y1(q2p)jx10jsinlp7r  * 


[qj* 


r  q,  ~  ~  f  ~  •) 

]  2  J1^2e)  X3  +  [^iCV^Nr03  pir 


;  2 

q-j  * 


|  -f  JL(q2?)  X9  +  |2Y1(q2P).]Xl0jsln  pir 
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APPENDIX  B 


REPRESENTATION  OF  THE  FOURIER  COEFFICIENTS 
ASSOCIATED  WITH  THE  INCIDENT  WAVE 


This  appendix  contains  the  explicit  representations  of  the 
Fourier  coefficients  that  define  both  the  shape  of  the  Incident 
traveling  waves  and  their  Fourier  expansion  in  the  circumferen¬ 
tial  direction. 


The  coefficients  a  define  the  time  history  of  the  incident 

P 

pulses  and  (see  Fig.  2  and  Sections  III.l  and  IV. 1), 

and  are  determined  from  the  usual  expression  for  the  Fourier  con¬ 
stants  .  They  are 


1 

T 


xx 


sin 


m  t 
P 


dt 


in  the  case  of  an  incident  dilatatlonal  pulse,  and 


T  sin  co  t  dt 
xy  P 


(B.l) 


(B>2) 


in  the  case  of  an  incident  distortional  pulse. 

Since  a„  defines  the  shape  of  the  pulse,  identical  values  of  a_ 
P  p 

will  be  obtained  in  both  cases  if  the  shapes  of  the  incoming 
pulses  are  identical. 


Referring  to  Fig.  2a  for  the  time  histories  and  associated 
parameters  of  the  rectangular,  triangular,  and  linear  rise- 
exponential  decay  pulses,  we  obtain  after  integrating  (B.l)  or 
(B.2): 
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where  cm  is  the  speed  of  propagation  of  dilatational  waves  when 
(B.l)  is  being  calculated;  and  is  the  speed  of  propagation 

of  shear  waves  when  (B.2)  is  being  calculated;  op  is  the  Lanczos 
factor  discussed  previously  in  Section  VI. 2. 

In  this  case,  any  two  of  the  following  parameters  may  be 

qA  O  T  C  t 

selected  arbitrarily;  (2^  )»  (zb")*  (~b)  *  remaining  one 

is  determined  by  the  relation 
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where  the  parameters  a  and  will.  In  general,  be 

prescribed  input;  the  parameter  k  is  selected  so  that  in  the 

♦■■t  me  j»i  1  CO?!f  AfTAAfl  Will  hflVA  dftCAVAd 

to  a  small  value  (i.e.,  e"k  «  1);  and  only  asm,  of  the  two 
c  T  c  t0 

parameters  may  be  chosen  arbitrarily,  the  other  being 

determined  by  the  relation 


(B.8) 


The  coefficients  defining  the  circumferential  expansion, 
in  Fourier  series,  of  the  incident  pulses  will  be  discussed  now. 
As  indicated  in  Sections  I1X.1,  III. 6,  and  IV.2,  it  is  possible 
to  obtain  representations  of  these  coefficients  in  terms  of 
Bessel  functions  of  a  complex  argument. 

An  integral  representation  of  the  Bessel  function  of  the 
first  kind  with  complex  argument  is  given  by  kef.  22  as 


(i)nirJn(z)  ■  el*  008  ^  cos  nQdfl  ,  (n  ■  0;  1,  2,  •••)  •  (B.9) 

'’o  '' 

Employing  the  trigonometric  identities 

cos  r©  cos  s©  ■  2  [cos(r  +  s)©  +  cos(r  -  s)©] 
sin  r©  sin  8©  ■  -  j  [c°8(r  +  s)©  ■  cos(r  -  s)©] 


in  conjunction  with  (B.9),  we  are  led  to  the  following  representa¬ 
tion  for  the  Fourier  coefficients  defined  by  Eqs.  (3.6.4),  (3.6.7): 

V  -  ■  (1>n+1 6  l,3{<i+«)Jn<i3>  -  (V)  [w^ + ]} 

(B.lOa) 
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ApO  “  e  q3[<1+«)J0(q3)  -  (l-€)J2(q3)]  (B.lOb) 

R  m  A  3/  f  t  /  _  \  ..  t  /_.  \  1  /•>  <A  < 

-pn  “  V  2  ;riri-2v'13/  Un-2VH3' J  \o.^j 

Cpn  “  ^  e  (~?b^)  [Jn+1^3^  "  Jn-l*q3*]  (B.iOd) 

CP0  “  (i)  e  <l3(“f,b1)Jl<q3)  (B .  lOe) 

P 

“pn  -W"*1  8  ^("^[WV  +  Jn-l(V  ]  <*•«*> 


The  preceding  coefficients  ere  of  course  complex*  For  use 
on  the  right-hand  aide  of  the  equivalent  real  matrix  equations 
discussed  in  Appendix  A,  ve  separate  (B«10)  into  real  and 
imaginary  parts  (denoted  by  superscripts  &  and  I,  respec¬ 
tively).  Thus,  recalling  the  definitions  of  the  complex  quan¬ 
tities,  q3  and  <  [Eqs.  (3.4.16)  and  (3.4.18)  ],  and  with  the 

other  barred  quantities  defined  subsequently,  ve  obtain 
Apn  ■  •  ^  i^n  «•  .in  53J 

AJn  ‘  *  3  [Apn  008  ’3  ■  8l“  %] 

Bpn  "  0  3  [®£n  «8  ^3  +  “pn 

»Jn  •  *  3  PpV  “•  <3  -  5pn  Ij)  (#  l 

“  e  "3  (“pn  cos  ’3  +  5pn  8ln  ’3] 

Cpn  *  6  3  [=J„  “•  *3  ’  S3] 

*V» 

“pn  -  8  3  [“£„  “»  S3  +  “pn  8ln  ’3] 

“pn  *  8  ^  (“pn  «•  S3  *  “pn  *in  «,]  - 
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where 


V  "  r[e(*n  ‘  ^2)  +  (l+e>1„  +  u -«>>,] 

Apn  "  iK1,,  '  &)  '  U+*Hn  -  a-«)I'2] 

*P°  ’  2{‘[J0(q3)  +  J2<‘»3>  j  +  a+OJ^V  -  (l-«)j](q3)J 

AP°  ■  +  J2'v]  -  a+ojJ(q3)  +  a-7)jf(,3)] 

•U- ♦  a-^J 
*  *  ^(vSx ♦  v:J 

^  -  V?<v] 

cp0  ■  £b[  V?<i3>  +  Vl<V  j 

♦  v„M 

5p»  ‘  ‘  i  yfoC  -  ¥’1]  • 
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where 


& m  .  i.  i^aL.\ 

\ca  b/  pir  V  2b  1  *  S 

P 

and,  when  n  Is  an  even  integer, 

\  •  <-l>“/2  ' 

^  -  (-l)n/Z  r[jnl(q3)j 

?1  -  • 

-  (-D^2 

IJl  «  f(,in^3>  +  Jn-1^  3 

(8.13) 

^  -  -  X-l)^2  f^i<l3> 

*[42(,3)  + 

&  *  -  (-1>R/2  +  Jn-2<’3>  ] 

I‘z  -  C-1>”/2  ?[jJ+2(v  -  j5.2(V  ] 

*«2  ■  ^1)n/2  fl^V  -  Jn-2«3>  ]  * 
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whereas,  when  n  Is  an  odd  integer. 


T  m  -  ✓  — 1  \  ^  _f  •!/-  *  1 

Xn  '  l'  TWJ 


^  -  <-l> 2  j 

^  ?{4h«3>  -  ] 


n+l 

(-D  2, 

* 

4  -  (-l>  2  2 


Jn+l(V  "  Jn-l(q35 


B±1 


^  1  *  ^  2  2  + 


stL 


£l  -  -  <-«  *  +  JU 


oil 


^  -  (-1)  *  ♦  4j(q,)] 


B& 


**2  *  <~l>  2  ?l,£+2<q3)  +  ^2<q3>] 


9±i 


*n2  *  (-1)  2  2  Jn+2(q3)  ~  Jn-2*q3*. 


B±I 


*n2  "  *  2  2[Jwf2*q3*  "  ^-2**3* 
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The  coefficients  defining  the  circumferential  expansion,  in 
Fourier  series,  of  the  stresses  and  displacements  associated  with 
the  incident  shear  wave,  are  given  by  Eqs.  (4.2.16)  and  (4.2.19). 
In  a  manner  similar  to  that  employed  in  the  case  of  the  incident 


- 1  •.«. —  *.L 

u  bawuai.  avitaj.  'fu  vv|  vtiwtjw  ww^4.*v4VMtow  m«.  v  wuwnti 


V  •»  w/Nnvn  r>  nM  nkl  m 


in  terms  of  Bessel  functions  of  a  complex  argument.  Employing 
Eq.  (B.9),  and  the  trigonometric  Identities  written  immediately 
below  it,  we  are  led  to  the  following  representation  for  the  co¬ 
efficients  defined  by  Eqs.  (4.2.16)  and  (4.2.19): 


A  -  -  U)n+1  e’1’4 
pn 


Jn+2(q4)  “  Jn-2(V 


Bpn  “  ^)n+1  e’i<l4 


Jn+2(V  +  Jn-2(V 


Bp0  “  1  e 


:<V  ] 


pn 


-  Qf)\  e  1,4  +  WV 

P 

~i(l. 


Dpn  *  (oH?)q4  e  4  Jn+l(V  "  Jn-l^V 


°P0  -  <$)  ’A  e  [J1<V 


(B. 15a) 

(B.15b) 


(B.15c) 


(B.15d) 


(B .  15e) 


(B.15f) 


where  it  is  noted  that  the  above  coefficients,  A  ,  B  ,  ...,  D  A, 

'  pn'  pn  p(r 

associated  with  an  incident  wave  of  distortion  are  defined  differ¬ 
ently  from  the  corresponding  coefficients,  Eq.  (B.ll),  associated 
with  an  incident  wave  of  dilatation. 


As  in  the  case  of  the  coefficients  corresponding  to  the  inci¬ 
dent  shear  wave,  we  separate  the  above  complex  coefficients  into 
real  and  imaginary  parts  denoted  by  superscripts  R  and  I,  re¬ 
spectively.  Recalling  the  definition  of  [Eq.  (3.4.16b)], 

and  with  the  other  barred  quantities  defined  subsequently,  we  ob¬ 
tain  expressions  similar  to  those  given  by  Eq.  (B.ll): 
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R  — o  —  — 

Av_  ■  e  A*  cos  q,  +  A 
pn  pn  n4  pn 


L  sl“  ?4] 
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,R  /*4  frl 
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*V 

"1a  r-i 
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<*s# 
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I  f— T  —  —  o  — 
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pn  pn  ^4  pn 


44  pn 


(B.16) 


III 
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with,  in  this  case. 


where  the  barred  and  tilde  quantities,  Xn2'  rnl'  are 

tions  defined  by  Eqs.  (B.13)  and  (B.14),  with  argument  q^ 
placed  by  q^. 


(B.17) 


func- 

ra- 
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